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ABSTRACT 



In the present work we review the most prominent lepton flavor violating 
processes (/i 67, /i 3e, (/i, e) conversion, M — M oscillations etc), in 
the context of unified gauge theories. Many currently fashionable extensions 
of the standard model are considered, such as: i) extensions of the fermion 
sector (right-handed neutrino); ii) minimal extensions involving additional 
Higgs scalars (more than one isodoublets, singly and doubly charged isos- 
inglets, isotriplets with doubly charged members etc.); in) supersymmetric 
or superstring inspired unified models emphasizing the implications of the 
renormalization group equations in the leptonic sector. Special attention is 
given to the experimentaly most interesting (/x — e) conversion in the pres- 
ence of nuclei. The relevant nuclear aspects of the amplitudes are discussed 
in a number of fashionable nuclear models. The main features of the relevant 
experiments are also discussed, and detailed predictions of the above models 
are compared to the present experimental limits. 
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1. INTRODUCTION 



All currently known experimental data are consistent with the standard 
model of weak and electromagnetic interactions (SM). Within the framework 
of the SM, baryon and lepton quantum numbers are seperately conserved. In 
fact one can associate an additive lepton flavor quantum number with each 
lepton generation which appears to be conserved. There are thus three such 
conserved quantum numbers Lg, and Lr each one associated with the 
lepton generations (e~, z/g), z/^), (r~, v^), with their antiparticles having 
opposite lepton flavor. It is in fact these quantum numbers which distinguish 
between the three neutrino species if they are massless. 

Most theorists, however, view the SM not as the ultimate theory of nature but 
as a successful low energy approximation. In possible extensions of the SM it 
is legitimate to ask whether lepton flavor conservation still holds. In fact in 
such gauge models (Grand Unifled Theories, Supersymmetric extensions of 
the SM, Superstring inspired models) such quantum numbers are associated 
with global (non local) symmetries and their conservation must be broken at 
some level. 

Motivated in part by this belief the search for lepton flavor violation, which 
began almost half a century ago (Hincks and Pontecorvo, 1948 [|T|, Lagarigue 
and Peyrou, 1952 0, Lokanathan and Steinberger, 1955 0, see also Frankel, 
1975 [Q) has been revived in recent years and is expected to continue in the 
near future. In the meantime the number of possible reactions for testing 
lepton flavor has been increased. The most prominent such reactions are 

/X ^ 67 (1) 
r — >■ 67 and r — >■ /i7 (2) 

/i — *• 66"'"6~, (3) 

r — »• ee^e~ , r — >• fie^e^ (4) 
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r — i> e/x"*"/! , r /i/i^/i (5) 
/w^e^, ^ 7r+/ie (6) 

(/i^e^) ^ (/U^e^) muonium — antimuonium oscillations (7) 

fi~{A, Z) ^ e~{A, Z) {muon — electron conversion) (8) 
Finally one could have both lepton and lepton flavor violating processes like 

(A, Z) {A, Z±2) + e^e^ {pPou - decay) (9) 

fi~{A, Z) — i> e^{A, Z — 2) {muon — positron conversion) (10) 

From an experimental point of view the most interesting reactions are (|l|), 
(0)) (H)) (ED (0)- ^^^^ report we will only briefly be concerned about 
the last two reactions. 

The problem of lepton flavor non-conservation is connected with the family 
mixing in the lepton sector. Almost in all models the above process can 
proceed at the one loop level via the neutrino mixing. However, due to 
the GIM mechanism in the leptonic sector, the amplitude vanishes in the 
limit in which the neutrinos are massless. In some special cases the GIM 
mechanism may not be completely operative even if one considers the part 
of the amplitude which is independent of the neutrino mass (Langacker and 
London, 1988 §, Valle, 1991 0, Gonzalez- Garcia and Valle, 1992 §). Even 
then, however, the process is suppressed if the neutrinos are degenerate. It 
should be mentioned that processes (|])-(|D cannot distinguish between Dirac 



and Majorana neutrinos. Processes (§) and ([T0|) can proceed only if the 
neutrinos are Majorana particles. 

In more elaborate models one may encounter additional mechanisms for lep- 
ton flavor violation. In Grand Unified Theories (GUT's) one may have addi- 
tional Higgs scalars which can serve as intermediate particles at the one or 
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two loop level leading to processes (|T])-(||). In supesymmetric extentions of 
the standard model one may encounter as intermediate particles the super- 
partners of the above. Lepton flavor violation can also occur in composite 
models, e.g. technicolor In fact, such models have already been ruled 
out by the present experimental bounds (see next section). 



The observation of any of the above processes (eqs. (P)- ([TO|) ) will deflnitely 
signal new physics beyond the standard model. It will severely restrict most 
models. It may take, however, even then much more experimental effort to 
unravel specific mechanisms responsible for lepton flavor violation or flx the 
parameters of the models. The question of lepton flavor non-conservation 
has been the subject of several review papers (Scheck, 1978 0, Costa and 

Vergados, 1986 ]12|, Melese, 
, Schaaf, 1993 HI). In the 



Zwirner, 1986 |Tg, Engfer and Walter, 1986 
1989 \K 



Heusch, 1990 fg, Herczeg, 1992 p 



present review we will focus our attention on recent theoretical developments 
of the subject. We will only cover the essential points of the experimental 
situation since we do not intend to dublicate the recent experimental review 
which appeared in this journal (Shaaf, 1993 |]16[)- Furthermore, the reader 
can find an interesting account of the early experiments by Di Leila il^, [18 



2. LEPTON FLAVOR VIOLATING PROCESSES 

We have seen in the previous section that lepton fiavor violation, if it occurs, 
can be demonstrated by many reactions (see eqs. (0)-(D). In this section 
we are going to examine the most basic features of the experimentally most 
important processes. 



2.1. The /i — 67 process 

As we have already mentioned this is the oldest and perhaps the best studied 
process. It was expected to proceed quite fast since the muon and electron, 
with the exception of their mass, are identical and possess identical electro- 
magnetic and weak interactions (for a historical review see Di Leila, 1993 
[O, Vergados, 1986 |]12|)- such early estimates the branching ratio was 



(Feinberg, 1958 ^) 
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which was puzzhng since it was aheady an order of magnitude larger than 
the experimental limit at that time. We will see, however, in sect. 3.1.1 
that we have additional suppression due to the leptonic analogue of the GIM 
mechanism. In purely left-handed theories the branching ratio takes the form 
(Vergados, 1986 pl) 

R=^[€'+4'f (12) 

with 

e^ = E^if^f^;r^' m,«m^ (13) 

Tfl-fXT 



(L) 

Vn 



(12)..*(12) 



m 



w 



Ml 



aln — ir- + b 



m 



w , 



Mj ^ niw 



(14) 



where U, 



gj^^ {u'^^P^ is the amplitude for producing the light eigenstate with 



mass rrij -C mw in the weak interaction of e" (/U~). f/*-^^^ provide the corre- 
sponding amplitudes for the heavy neutrino components. For neutrinos less 
than IMeV the term r/^^^ is negligible. The term r/j^'' also becomes negligible 
for very heavy neutrinos. 

It is clear, therefore, that lepton flavor violating processes are suppressed. 
Precisely how suppressed one does not know. Thus, experimentalists have not 
been detered from pursuing such hard experiments. For purely experimental 
reasons only positive muons have been considered since, among other things, 
they do not undergo capture by a nucleus as negative muons do (see sect. 
2.2). The experiment consists in the simultaneous detection of a photon and 
a positron moving essentially back to back with momentum p ~ m^c/2. 

In addition to accidental ej events, which can be minimized by shielding, 
the main source of background is radiative muon decay i.e. 
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/i+ e+i/eZ^^7 (15) 

in the kinematic regime in which the neutrinos carry away very httle energy. 
In spite of the heroic experimental efforts, this process has not been observed. 
The best experimental limit 

R = ^ e+7) ^ ^ g ^ ^ 

r(/i+ e+z/gZ/^) 



has been set by LAMPF (Bolton et ai, 1988 [Q) using the crystal box de- 
tector. This is almost an order of magnitude improvement over the previous 
record, R < 1.7 x 10^^'^ (Kinnison et ai, 1982 |21|). The proposed limit by 



the new LAMPF detector MEGA (magnetic spectrometer with large solid 
angle) is 

R ~ IQ-^^ (17) 

the final sensitivity depending on the availability of LAMPF beam time 
(MEGA collaboration |||, Cooper et ai, 1985 [||). 



In the same experiment a limit on the branching ratio for two photon emission 
was indirectly set 

R < 7.2x IQ-^^ (90% CL) (18) 

For comparison we mention the corresponding limits for r decay obtained by 
the CLEO collaboration (Bean et ai, 1993 |4|) 



R = ~^ < 4.2 X 10-*^ (90% CL) (19) 
(r all) 

In the intermediate neutrino mechanism the branching ratio is still given by 
eq. ([T2|) with the modification Uej — > Urj provided of course that we ignore 
all chanells other than /i^z/^z/,- in the denominator (otherwise the r- 

mass dependence is complicated). Even in this case, however, the branching 
ratio is expected to be larger crudely speaking by a factor {m^^ /m^^Y due 
to the large coupling of the r to the heaviest neutrino mass eigenstate. The 
branching ratio for r — > 67 is i? < 1.2 x 10~^ obtained by Argus (Albercht 
et ai, 1992 |^). For a complete list see Depommier and Leroy, 1993 p6| . 
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2.2. The /i — > eee decay 

In principle, every model, which allows /i ^ 67 to proceed, will also allow 
fi — > ee+e~, the only difference being that now the photon can be virtual de- 
caying to an e~^e~ pair. In such models one expects /i — > 3e to be suppressed 
by a power of a i.e. 

i4 - — - 10-' (20) 
i?(e7) 247r ^ ^ 

One can construct models, of course, in which the opposite is true (see 
Bilenky and Petcov, 1987 P3)- Furthermore, since the photon is now vir- 



tual, one may have a contribution from the EO and MO form factors [|T2 
which vanish at = 0, i.e. for real photons. We should emphasize, however, 
that there exist models which allow fi —>■ 3e but forbit fi —>■ 67. Such are e.g. 
models in which lepton flavor can be mediated by Higgs scalars which are 



doubly charged |T2[. In fact, in such cases fi —>■ 3e can proceed even at the 
tree level. In short, if lepton flavor violating mechanisms exist, /i — 3e has 
a better chance of being allowed. 

In addition to the above theoretical considerations, /i — 3e offers a number 
of experimental advantages as well. One can take advantage of the three 
charged particles appearing in the final state to reduce the background by 
a variety of timing and kinemetic constraints. The chief required qualities 
for the detector are: good energy resolution, time resolution and precise 
vertex construction. A detector which meets well these specifications is the 
SINDRUM at PSI [^. The best upper limit obtained by such a detector is 

R = r(/i+ ""^""'p < 1.0 X 10-12 (90% CL) (21) 
It is worth noting that, this limit is a two order of magnitude improvement 



over the limit obtained by the same group (Bertl et al, 1985 ^9]), which in 



turn was an order of magnitude improvement over that of the Dubna group 



(Korechenko et al, 1976 ||3y]), which had stood for about 10 years. It is 
obvious that, such experiments should be encouraged to continue. At least 
they will provide supplemental information to — > 67. 
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2.3. The (yU, e) conversion in the presence of nuclei 



When a negative muon stops in matter, it finally forms a muonic atom cap- 
tured by a nucleus with a radius of about 200 times smaller than that of the 
usual atom and a binding energy in the KeV regime for light and medium 
nuclei (for heavy nuclei the binding energy is of the order of a few MeV see 
below sect. 4.2.1). After it cascades down to the Is orbit by emitting x-rays, 
it eventually disappears by decay in fiight 

fi~ e"z/eZ/^ (22) 

or by capture by the nucleus 

^I-{A,Z)^u,{A,Z-l) (23) 

The former dominates for light nuclei but it is negligible for heavy nuclei 
[Z > 10). If lepton fiavor is not conserved, one can encounter the process 

fj^~{A, Z) — s> e~ {A, Z) {muon — electron conversion) (24) 

The (yU~, e" ) conversion in the presence of nuclei is the most interesting lepton 
fiavor violating process from an experimental point of view. The reason is 
that, the detection of one particle is sufiicient. No coincidence is needed. The 
energy region of the produced electron is almost background free. There are, 
however, some sources of backgound which have to be taken into account. 
The first is the muon disintegration in orbit mentioned above eq. (^2]). The 
electron energy in free muon decay is below 53MeV". The bound muon decay, 
unfortunately, yields electron in the high energy tail. These electrons can be 
mistaken with the interesting (/x, e) conversion produced electrons due to the 
finite energy resolution. Good energy resolution is thus a critical and crucial 
feature of a detector. 

The second source of backgound is radiative muon capture 

^i-{A,Z)^u,^{A,Z-l) (25) 

followed by pair production. This source of background can be eliminated if 
one considers only (/i, e) conversion leading to the ground state of the final 
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nucleus. Then, by a judicious choice of the target nucleus the maximum 
electron energy following pair production in ( P^D can be much smaller than 
that of the monochromatic electrons of (/i, e) conversion |^. The backgound 
free region can be as big as 2.5MeV. 

The third source of background, which can be minimized by shielding, comes 
from electrons produced by cosmic rays. Fortunately, this kind of backgound 
can be studied experimentally, since it is easy to accumulate good high statis- 
tics data on cosmic rays during the beam off periods. 

Another troublesome background comes from radiative pion capture which 
is followed by pair production. This can lead to copious electrons in the 
interesting energy region. Fortunately, the life times for muonic atoms are 
quite long (70ns or longer). So, conversion events can be seperated from 
the prompt background. It is, however, important to eliminate most pion 
contamination in the beam with electromagnetic seperators. 

It is clear from the above discussion that, the detectors must have good 
energy resolution and large solid angle. Such detectors with good detection 
efficiency must meet conflicting requirements. One usually tries to strike a 
compromise between the expected sensitivity, which rises with Z, and the 
signal to background ratio which drops with Z. It is also important to ensure 
that the coherent production exhusts a large fraction of all the conversion 
electrons. Fortunately this happens to be the case (see sect. 5.2.2). 

The flrst most efficient such detector for (/x, e) conversion is the Time Pro- 
jection Counter (TPC) developed at TRIUMF (Bryman et ai, 1993 |2[). 
Using a target of Ti, dictated by the above requirements as well as addi- 
tional experimental constraints, the TRIUMF collaboration has obtained 

(Ahmad et ai, 1988 ||33[). Using a Pb target the same group has obtained 

= < 4.9 X 10"^° (90% CL) (27) 
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At the same time the experiment found no positron candidate for process 
(|TU|) with excitation energy in the final nucleus below a few MeV. This has 
led to the following limits 



R = T{ij,-Ti e+Ca{gs)) < 0.9 x 10"^^ (90% CL) (28) 

The total branching ratio depends, of course, on the final nucleus excitation 
energy. For a giant resonance distribution with a mean value of 20MeV^ the 
authors deduced the limit 



R = T{fi~Ti e+Ca{E < 2meV)) < 1.7 x 10"^° (90% CL) (29) 

Another detector with high sensitivity is the SINDRUM II Spectrometer 
(Badertscher et al, 1991 @]). With this detector, during the test run, a 
marginal improvement was achieved for the (/i, e) conversion branching ratio 

= < 4.4 X 10-^^ (90% CL) (30) 
while the (/i~, e"*") has improved by a factor of 2, i.e. 

R{fi-Ti ^ e+Ca{gs)) < 5.5 x 10"^^ (99% CL) (31) 
The next goal set by SINDRUM g is 

R ^ 3x 10"^^ (32) 



and by MELC (Djilkibaev and Lobashev, 1992 §^) is 



R ^ 10~^^ (33) 

It is clear from the above discussion that, {fi, e) conversion has definite ex- 
perimental advantages. The proposed limits may meet the predictions of 
realistic models. Furthermore, (/z, e) conversion, like /i — 3e discussed in 
the previous section, may occur in a number of models which do not lead to 
{fj, 67) such as e.g. those involving the box diagrams (see sect. 3). Due to 
the fact that, the various mechanisms (discussed in sect. 3) lead to different 
A and Z dependence, (/x, e) conversion, if it is ever observed, may be able to 
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shed light even on the detailed mechanisms for lepton flavor violation. For 
this reason it will be discussed in detail in sects. 4 and 5. 



2.4 Muonium-antimuonium oscillations 

The muonium atom M = (yU+e~) has interesting electromagnetic properties 
compared to positronium. In the presence of lepton flavor changing inter- 
actions transition between muonium and antimuonium, M = {fi~e~^), can 
occur which are the analogue of the well known oscillations in the — 
system. The only difference is that the expected oscillation period is much 
longer compared to the life-time of Additional complications occur due 
to the fact that the degeneracy of M and M may be destroyed in an external 
magnetic field (> O.OSG) or due to interactions with matter. 

The amplitude for M — M oscillations takes the value of a typical four- 
fermion weak interaction multiplied by a scale factor (Vergados, 1986 



[0), which depends on the gauge model, and represents the lepton flavor 
violating parameter. It is this parameter which we expect to extact from 
the oscillation measurements. Even though the oscillation time is inversely 
proportional to n^, 

r ~ 1.3 X lO'^n-^ (34) 

in the actual experiment one measures the probability for M decay which 
takes the form 

Pm ~ 2.6 X IQ-^n-^ (35) 

Further reduction can occur in the presence of a magnetic field, as we have 
already mentioned above. 

An important step towards experimental detection of M — M oscillation was 
achieved after the development of techniques to produce thermal muonium 



in vacuum (Marshall et al, 1982, 1988 and Huber et al., 1990 [36|). Another 
important step was the ability to measure a low energy atomic positron in 
coincidence with the electron of /i~ decay (Mudinger et al, 1988 and 
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the capacity of SINDRUM I (Jungmann et ai, 1989 [^) spectrometer to 
increase the electron sohd angle by roughly a factor of 300. 



In spite of the above important steps, the sensitivity of M — M oscillation 
experiments is limited not only by the smallness of rix, but by the additional 
factor of 2.6 x 10~^ of eq. ( ^5|) (Schaaf, 1991 |]l6l). It may, however, become 
of interest in some special models, especially those for which the process can 



occur at tree level (Herczeg and Mohapatra, 1992 |3^, Vergados, 1986 |T2| 



The highest sensitivity has thus far been achieved by LAMPF (Matthias et 
ai, 1991 [0). No candidate events were found in the experiment which led 



to the upper limit for 



Pm < 6.5 X 10" 



(90% CL) 



(36) 



This leads to 



Ur < 0.16 



(37) 



The SINDRUM collaboration plans to reach a sensitivity of 10 by 1993 
with the ultimate of 10~^^ one year later. 



2.5. Lepton flavor violating meson decays 

At first sight, the best such process seems to be 7r° — > /x^e^. It turns out, 
however, that the branching ratio for this reaction is small, since it has to 
compete not against a weak competitor but against the electromagnetic decay 
7r° — i> 27. Thus, the most prominent such decays are 

Kl ^ /i^e^ (38) 

K+ 7r+/ie (39) 

in spite of the fact that the experiments have to be done in flight with not 
so intense beams. Such experimental efforts are expected to intensify, if 
any of the planned meson factories are ever constructed (TRIUMF, Euro- 
pean, Moscow). Historically, such processes have been favored, because of 
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some kind of prejudice in favor of generation number conservation. In other 
words, a change of generation in the leptonic sector may be compensated by 
a corresponding change of generation in the quark sector. 

For reasons analogous to those mentioned in the discussion of fi —>■ 3e, the 
reaction (^) appears to have some advantages. In other words, it has the 
advantage of abundant particle identification, which may be used to discrim- 
inate against background i.e. vr+vr" pairs being mistaken as fi~^n~ pairs. The 



best limit ever set for the branching ratio of process (39) is 



R={K+ ^ 7r+fie)/{K+ all) < 2.1 x 10~^° (90% CL) (40) 



by Lee et ai, 1990 This is almost an order of magnitude improvement 
over the previous limit < 1.1 x 10~^ (BNL E777 experiment, Campagnari 
et ai, 1988 [0). The branch /i+e" rather than the branch fi~e^ was selected 
partly due to the generation argument mentioned above (s ^ fi^,d ^ e~) 
but mainly due to the fact that the positron background is more formidable 
than the electron background. 

At present, the most sensitive limit comes for the reaction (^8|), i.e. 

R = < 3.3 X 10-11 (90% CL) (41) 

— ^ all) 



which has been set by Arisaka et ai, 1993 Once again, this is an order 



of magnitude improvement over the previous experimental limit, 

R < 3x 10-1° (42) 
which has been set by the BNL E791 experiment (Cousins et ai, 1988 



3. LEPTON FLAVOR VIOLATION IN GAUGE THEORIES 

We have pointed out in the introduction that, in the standard model of elec- 
troweak and strong interactions neutrinos remain strictly massless and lepton 
flavor is automatically conserved as a global symmetry of the Lagrangian. In 
most of the extensions of the standard model, however, there are various 
sources of lepton flavor mixing and processes like cy, /i — > eee, /i — > e 
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conversion etc., which occur at the one-loop level. In particular, lepton flavor 
non-conservation appears in the cases where these extensions predict non-zero 
neutrino masses or an extended higgs sector. Natural candidate theories are 
Grand Unified Models HH] and Supersymmetric Theories HH]. In this sec- 



tion, we are going to summarize briefiy the most important extensions of the 
standard model. We start with a general discussion of the neutrino mass 
mechanism when a right handed neutrino is included. We further discuss 
models with additional scalar particles (singly charged isosinglets, doubly 
charged isosinglets and isotriplets) , and finally we present a brief overview 
of flavor violating effects in Supersymmetric and String motivated Grand 
Unifled models. 



3.1. Minimal extensions of the standard model 

3.1.1. The right handed neutrino. The most obvious way to extend the 
standard model is to include the right handed neutrinos. When the right 
handed neutrino is present, a non-zero Dirac mass term is possible in the 
theory and a "Kobayashi-Maskawa" leptonic mixing matrix appears, which 
gives rise to flavor violations. 

Moreover, since neutrinos are electrically neutral, they can in principle have 
Majorana masses, violating the total lepton number by two units. In this 
latter case, new processes may also occur, namely the neutrinoless double 
beta decay (/9/9oi/-decay) and muon to positron conversion in the presence of 
nuclei /i" + [A, Z) e+ + {A, Z -2). 

Thus, in the general case one may have the following new mass terms in the 
Yukawa sector of the theory 

Cmass = I'lmvl'^ + Nl'm^Nl + vlmj,Nl + Nl'mlv'j^ (43) 

where, vl = (z/°, z/°, 1/°)^ and = (iV°, iV°, iV°)K, are the left and right 
handed neutrino weak eigenstates. With 

'^'n=C{ulf , Nl'^ = C{N'^f (44) 
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we denote the conjugate fields, while 171^,1710, iTiJy,mN are 3 ® 3 matrices. 
Thus, the most general neutrino mass matrix is the 6 6 matrix 



(45) 



The matrix (^) can be diagonalized by separate left and right unitary trans- 
formations. Assuming that the masses of the left handed neutrinos are much 
lighter than those of the right handed ones, and labeling their eigenstates 
i/jx, and iVjij, respectively, the transformation for the left - handed fields is 



For the right-handed fields the transformation is 




(46) 



(47) 



In eqs. ( ^Bj ) - (^7D, = Aj(z/) and Ajy = Aj{N) are diagonal matrices of 
arbitrary phases, while a = aj and (p = (fj are phases related to the CP- 
eigenvalues of the neutrino eigenstates as 

{CV)iyj{CV)-^ = e-^'^Vj , {CV)Nj{CV)-^ = e'^'^Nj (48) 
Due to the presence of neutrino mass, the charged leptonic currents remain 



no longer diagonal. Thus the left handed current becomes |T2|, 



= -2e°7^z/° + h.c. = -2{eLl^U^^VL + clI^U^^Nl) + h.c. (49) 

with 

U^^ = e-'^^{Sl)^Sfe-'^^ (51) 

where is the charged lepton mixing matrix. The charged current involving 
the conjugate fields becomes 
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if^y = -2[unl,e-iU'Ye% + NRl,e^nU'Ye%] (52) 

Thus, the phases aj, ipj are in principle measurable. They appear in processes 
where the conjugate fields are present {jSPoi, - decay and fj,~ — >• e"*" conversion). 
The right-handed current is modified analogously 



= 2{eRj^U^'un + cr^^U'^Nr + h.c.) (53) 



where 



^21 ^ gi(A.+A.)(^e)T^21*g-ia (54^, 



f/22 = e*(^=+A«)(5i)tSf*e-'^ (55) 

Notice that, in the presence of right-handed currents the S"^^ and S*^^ matrices 
are involved. If the right-handed neutrinos are absent, in the above formulae 
jji2^jj2i g^j^^ ^22 g^j^g zero. U^^ is the leptonic analogue of the Kobayashi- 
Maskawa mixing matrix. 



The above considerations apply in most Grand Unified models (GUT's) P5| 
where non-zero neutrino masses arise naturally. In Supersymmetric GUT's in 
particular, motivated by the observed merging of the Standard Model gauge 
coupling constants, there has been a revived interest in determining the low 
energy parameters of the theory, including the unknown neutrino masses and 
mixing angles |^^, in terms of a few inputs at the GUT scale. The general 
strategy in these approaches is to use the minimal number of parameters at 
the GUT scale, so as to have the maximum number of predictions at mw 
Ultimately, one hopes that this minimal set of parameters at the GUT scale 
may be justified in terms of a more fundamental theory, such as the String 
Theory. The advandage of such a procedure is that there are many direct 
or indirect constraints on the neutrino masses from the rest of fermions. For 
example, the Dirac neutrino mass rriy^ is usually related to the up-quark 
masses in most of the GUT models. It thus appears challenging to utilize 
all possible such constraints in the neutrino mass matrix, in order to make 
definite predictions for the as yet elusive neutrinos, which will then be checked 
by experiment, this way supporting or excluding such GUT scenarios. In the 
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next sections we are going to use particular predictive models for fermion 
masses in our estimations of the flavor violating branching ratios. 



3.1.2. The extended higgs sector. In this subsection, we are going to 
review the basic features of minimal extensions of the standard model based 
on the introduction of new scalar particles [|^, 0, ^ which are consistent 



with the gauge symmetry. Neutrino masses are generated in this case without 
introducing any right handed neutrino. Possible scalars, which couple in a 
renormalizable way to leptons, are additional higgs doublet fields H^'''^ = 
(2, i), = 1, 2, a simply charged scalar singlet field S~ = (1, —1) ^ 
a doubly charged state A^"*" = (1, —2) as well as a triplet T = (3, —1). 
All the above states can couple to leptons and create diagrams leading to 
lepton flavor non-conservation. Their Yukawa couplings are |^ 

6C = \,ld'^RS + d,,IiJ-^^]R + ■ie.ijA* + h.c (56) 

It is possible to conserve lepton number in the above Yukawa Lagrangian, if 
in the above fields we assign the following lepton numbers: 

L{S) = 2, L(T) = 2, L(A) = 2 (57) 
But lepton number can be violated explicitly by cubic as well as quartic 



terms of the above fields 50 



/z(i/°*, -H~*) [^1)S-+ fi\H'*, -H-*)^ ( ^0 ) + - (58) 

where H' is a second doublet and /i, /i', .•• are the qubic coupling mass param- 
eters. (Notice that couples antisymmetrically to isospin 1/2 particles.) 

A generalization of the above model involves, one additional higgs doublet H', 
the additional isosinglet fields S""^ and the singlet scalars where a, b stand 
for e, fi, T. Here, rather than breaking the lepton flavor explicitly |4^, ^ 



we prefer to follow the approach of references [^2|, |5^ and consider models 
that above some scale V, have a global abelian lepton flavor symmetry G = 
f/(l)e X f/(l)^ X U{l)r. At the scale V, the group G is broken spontaneoulsy 
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Figure 1: One loop contribution to the neutrino masses in the model with 
an extended Higgs sector. 



by vacuum expectation values (vev's) of the singlet scalars giving rise to 
Goldstone bosons: ^""^^ = e*'^'=''/^<=f , ^"'^ = e^^<^^l^<^^ and ^^'^ = e^-^"^/^"^. Thus, 



the new interactions in the Lagrangian are p2l p3 



tI nabf 



(59) 



Neutrino masses are generated by the one-loop diagram shown in fig. 1. The 
loop calculation for this graph gives: 



mab 



Kb gba < ^ > ^ l-n 



167r2 



M| - ml 



w 



(60) 



In the above, ma, mb are the masses of the charged leptons, Ms is the mass 
of the heavy charged singlet 3°"^ which appears in the loop, while v is the 
vev of the standard higgs doublet. We have also assumed that, the second 
doublet circulating in the loop, which does not develop a vev, has a mass of 
order m^y. 

In the most general case, where all the singlet fields acquire vev's, the 



neutrino mass matrix can be parametrized as follows |53 



/ 



mo 



tanO 




(61) 



tanO 
y coscj) sincj) 

In the above matrix mo sets the mass scale and is given in terms of the 
various parameters entering eq. (BO) from the formula 



mo 



1 v{Vl + V^j'^ 

-m^/n 



167^2 M2 - ml 



w 



(62) 
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where we have made the approximations + ^ rrir and + mg ~ m^, 
while 

Vab = Xab~gbag'a<^"'> (63) 

Furthermore, since we know nothing about the couphngs g'^, in the last equa- 
tion we have used the approximation g'^ ^ g'^ ^ g'^- ^ g' ■ Finally, we have 
defined 

tancj) = ^ (64) 

''er 

and 

tanO = ^^^'^^ sincj) (65) 

In order to diagonalize the above matrix, guided by phenomenological rea- 
sons, we make the natural assumption that sin2(j) << 1. In this case, the 
eigenmasses are found to be 



rriy^ ^ — -mQsin26 sin2(f) (66) 

m^^ ~ %r + -mQsin29sin2(j) (67) 

cosO 4 



m 



^3 



rriQ , 1 



cosO 4 



H — mQsin26 sin2(j) (68) 



Furthermore, we get two different sets of eigenstates, depending upon whether 
sincp << 1 or coscj) << 1. Thus, in the case where coscj) << 1, the eigenstates 



are related to the weak states as follows 531: 



Uf, ~ cos9{sin(j) — cos'^9cos(f)sin2(j))h'i 

+ —j={sin(j)sin6 — cos0) z/2 + —j={sin(j)sin6 + cos0) z/3 

z/^ cos9{cos(f) — cos^0sin(j)sin2(j))ui 

+ —j={cos(j)sin6 — sin(f))h'2 + —j^{cos(j)sin6 + sm0)z/3 
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-sinOvi H — ■i=cos9{l 
v2 



+ 



1 

71 



sin2(f)sin6)i'2 
cos6{l + sinO sin2(f))v^ 



(69) 



A simpified version of the above model arises when one considers the par- 
ticular value coscj) = 0. This case corresponds to the particular choice 
< $er >= 0, which exibits an exact Le-^+r symmetry. In this latter case, 
one finds a zero mass for the first neutrino and two completely degenerate 
states for the other two, i.e. 



52, 53 



0, 



cosU 

The diagonalizing matrix is given by 



mo 



m 



1^3 



mo 
cosO 



U 



( cosO 


y-sin6 



sind 


sind 


\ 


V2 


V2 


1 


1 




'V2 


V2 




cosd 


cosO 


) 


v/2 


V2 



(70) 



In the last two equations there are only two parameters, namely, the mass 
scale mo and the angle 6. If we adopt the MSW-solution for the solar neu- 
trino problem, and a natural hierarchy between the mj,^ and rriy^ neutrino 
mass eigenstates, both parameters can be fixed uniquely. ( We mention here 
that an alternative solution to the solar neutrino "puzzle" is based on the 
assumption of a large neutrino magnetic moment; for a recent review see 
Pulido 0.) 



The model under consideration gives the following formulae for the various 
oscillation probabihties 

P{ue ^ Vr) ~ -sm^26'(l - -sm^vr-^) 

2 2 /23 



P{vr — * v^) ^ cos^Osin^Tx- — 



P{h'e ^li) ~ sin Osin n 



L 

23 

L 

^23 
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Here, the short oscillations have been averaged out. We notice that, the above 
oscillation probabilities are expressed solely in terms of two parameters, 9 and 



Z23. If the MSW-effect is interpreted through z/g 
the experimental data one finds that 



oscillations, then, from 



^ (1.79 - 3.46) x lO'W, sinO ^ (0.71 - 1.10) x 10^' (71) 
see also |^6|). We notice, however, that the above model can also accom- 



modate a relatively "large" neutrino mass without creating any particular 
problem in low energy phenomenology. As an example, we mention the solu- 



tion given by the proposed model 53 1 to the VJKeV neutrino "puzzle 



Other neutrino mass mechanisms have also been proposed in the context 
of Grand Unified Theories. Of particular interest is the Witten mechanism 
3^, which is possible |Q in all GUT- models which include the right hand 



neutrino in a larger representation. The advandage of this mechanism lies 
in the fact that the Higgses need not belong to large representations of the 
corresponding symmetry. Other mechanisms giving mass at the two loop 
level have also been discussed in models predicting large magnetic moments 



3.2. The flavor violating decays 

We are now in a position to discuss the various violating processes and obtain 
numerical expressions for the corresponding branching ratios. The processes 
under consideration are the flavor violating processes yU 37, /i 3e, 
(/i~, e~) conversion as well as the lepton number violating processes (/x~, e"*") 
and /5/3o!/-decay. 



3.2.1. The ^ 67 and /i — > 3e decay rates. We start with the familiar 
H —>■ e'j decay. From the theoretical point of view, this decay, due to its 
importance, has been discussed extensively in the literature the last two 
decades [Q. The essential features and experimental limits of this process 



have been discussed in sect. 2.1. This decay violates lepton flavor and can 
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proceed through diagrams which involve massive neutrinos or new scalars 
predicted in various extensions of the standard theory. 



The most general form for the on-shell (g^ = 0) amplitude for /i — 67 is 
given by 

M{fi^ e'j) = e {/ei + /A/i75)«em^o-p^g^e''/i (72) 

where e'^ is the photon polarization vector and a = |[7p, ju]- Once the form 
factors /ei, /mi are given, in a certain gauge model one can easily obtain 
the branching ratio with respect to — e + z/g + i/^ decay i.e. 



\fEl\'+\l 



- 67) = 2An\2na) , 'f^^' (73) 

In the mass mechanism there are basically four diagrams contributing to the 



amplitude [|12|. We distinguish the following cases: 



i) Left-handed currents only: 
The form factors in this case are 

\fEi\ = \fM.\ = ^^^,{v!:' + vk') (74) 

where f]^^, rj^^ are lepton violating parameters defined for light and heavy 
neutrinos, respectively, as 



= (75) 

'ii'W 

Vk' = Efif'cif ^^(.ta^ + t) (76) 



with a = 2 and b = —3. The above parameters in eq. ( |f6D include the mixing 
which is model dependent. In the case of a particular mass matrix ansatz for 
the SO (10) model, given for example in ref. |6l[], we have 



r/^^ ^ 5 X 10~^°, T]^^ ^ (77) 
which lead to a suppressed branching ratio 
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R ~ 10-^^ (78) 
Thus in left-handed theories — > 67 is unobservable. 

a) Right-handed currents (R-R couplings): 

In this case, we have further suppression due to the presence of the interme- 
diate boson Wr, which is assumed to be much havier than its left-handed 
partner 

f^V = A;<l (79) 

Assuming that the mixing C, of the two bosons Wl, Wr is small enough 
(C ~ 0.1/c), one obtains the following lepton violating parameters 



= k^Huf^uT'^ (80) 

if." = *^i:t'if'C*$('"»^+'') (81) 

which again lead to an unobservable branching ratio R < 10~^^ depending 
on the precise value of k. 

Hi) Left - Right mixing: 

This is the most flavorable case in the mass mechanism, since the obtained 
amplitude has different structure from the previous ones. We get 

rLR{RL) _r.GF ml ( LR{RL) , LR{RL)\ (^r,^ 
/ei,mi ~ " l^^ J ' 

where refer to the L(R)-handcd coupling for the iwW{euW) 

vertex and to the R(L)-handed coupling for the euW {fivW ) vertex. Their 
expressions are given by 



= CE^^f^*^ (83) 
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Figure 2: ^ e'y via the singlet field 5". 

V^' = CT.U^U^*- (84) 
# = CEt/if'c/r-;^ + (85) 

V^' = Cj:uT'uP*^(aln^ + b] (86) 

However, even in this most favourable case, in the mass mechanism the 
branching ratio turns out to be very small: 

From the above discussion it is clear that, the neutrino mass mechanism is 
in all cases negligibly small due to the tiny mixing and since rrij -C my/ or 

iv) extended higgs sector: 

It is possible to avoid the suppression mass mechanism by using an expanded 
Higgs sector. In sect. 3.2, we discussed such a possibility by introducing 
the isosinglet S~ . With the diagrams given in fig. 2 one can compute the 
branching ratio for the fi —>■ cj decay Q], which is given in terms of the 
mass Ms of the isosinglet. One finds |6^ 

which results to a bound for the singlet mass Ms 

Ms > MGeV X [10\X^Are)'^^] (89) 

• Another interesting process, which can occur with more mechanisms than 
/i —>■ 67, is the process n —>■ eee. Firstly, it can occur in all mechanisms which 
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allow /X — >• 67 decay, provided that the emitted photon is virtual. Secondly, 
this process can be mediated by the neutral boson Z, which finally decays 
to an 6"'"6~ pair. Finally, we can include now box diagrams |T^. Similar 
diagrams are also generated in the case of the Higgs singlet S. 

In the mass mechanism we can have contributions from all classes of dia- 
grams. We discuss here the most important ones. 

1). For purely left handed theories by comparing fi eee branching ratio 



with the corresponding one for /i 67 decay, one obtains [|12 



R ( i^Zl^] ~ 1.4 X 10-2 (90) 

The Z and W diagrams give a larger contribution than that of the photonic 
ones in the case of light neutrinos RzifJ' 36) ~ 9 x 10"^ R'y{^ 3e), 
RwifJ' 36) ~ 1.8 X 10'^-R7(;U 3e), but still far from the experimental 
limit. 

2). With the presence of the right handed currents, the only important 
contribution comes from the left - right mixing in the photonic diagrams. 
We get 



with 



where 



i? ~ 2.7 X 10"^ Ir^^^l^ (91) 



^LRiRL) ^ ^LRiRL) ^ ^LRiRL) ^gg^ 

Using the values of the previous sections for the neutrino masses and mixings 
one obtains 

R{l^ 36) < 10-3^ (94) 
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3). We can also introduce here the singlet S~ to avoid the mass mechanism 
suppression. Again three classes of diagrams are generated |T2|. The most 
important contribution comes from the photonic ones. The decay rate is 
given by 

r5(/^-3e) = ^mJ-|/Mi|^ (95) 
and the branching ratio by 

Rsifi 3e) ~ 0.73 x lO'^Rsifx 67) (96) 

3.2.2. (/i~ — e^) conversion decay rates. The next flavor changing pro- 
cess we are going to consider is the (/i~,e~) conversion in the presence of 
nuclei. The diagrams contributing to this process are similar to those we 



have considered for fi — > eee decay (see ref. |jT2|) 



Starting again with the photon diagrams, we write down the decay rate which 
is 



^ = 8rnl^a^^^eo\ME\^ (97) 

where = I /so + /miP + |/a/o + /eiP and \ME\ is the nuclear matrix 
element. (A detailed discussion will be presented in sects. 4 and 5.) For a 
wide range of nuclei and assuming only ground state transitions the matrix 
element \ME\ — > Mgg^gs lies in the range 

\ME\ Mg.^gs ~ (0.2 - 0.9)^2 (98) 

In the rest of this section we assume an average value \ME\ ^ 0.5Z^. The 
branching ratio is obtained with respect to ordinary muon capture fi~ + 
{A, Z) ^ Vfj^+ {A, Z — 1) decay rate which is ||63|, |6^ 



r(/i ^ u,) =ml — GlZtffZ [F^ + 3Fl + F| - 2F^Fp] /(A, Z) (99) 
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where Fy, Fa, Fp are form factors and the Primakoff 's function f{A, Z) for 
nuclei with A 2Z has the value f{A 2Z, Z) k, 0.16. Thus 

where is related to the Icpton violating parameters already defined in 
previous processes. It is obvious that, even in the L-R mixing where 
gets its maximum value ~ 9 x 10~^^, the contribution remains small. Thus, 
for light neutrinos 

RJlx- ^ e-) ~ l^-^^Z ^ (101) 

ml 

For the Z-diagrams the basic contribution arises in the L-L case. We thus 
get the branching ratio 



Rzin~ e-) ^ oo 2 V ^ ^7; V I'^zl ~ 0.34 x lO-^zT^ 



327r2 ml 8Z f{A,Z) ' ^' ' ^' 

(102) 



where 



■3 ,< > 



^1 ~ U + ^ + ~ 4 X 10-^3 (103) 

The quantity < > represents the effective quark mass. In the case of box 
diagrams, the L-L currents are also the most important and the branching 
ratio is found to be 

Rwil^- ^ e-) ^ 0.6a''Z^(ri'^ - Ar^^f (104) 

ml 

where 

= f 1 + /nl^) r;,-- + r;^^ (105) 

Here k (a) stand for the down (up) quarks. If we ignore the mixing in the 
quark sector and take ~ ma ~ 4:XlO~^GeV, we get Rw ~ ^O'^^ZEePe/ml 
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for light neutrinos, which is better that the previous cases. The presence of 
the right-handed currents in the box diagrams again appear to give neghgible 
contribution. 

Although fi — e conversion turns out to be negligible in the above context, 
as will be discussed in the following sections, one finds that this process is 
enhanced in super symmetric GUTs, when renormalization group corrections 
are taken into account. Let us, however, mention here that, in theories with 
extended fermion sector and a new neutral gauge boson Z', one may also 
have a significant impact in the /i — e conversion. It is argued that in 



such theories, the present experimental limits on the above process give a 
nuclear model independent bound on the Z — e — fi vertex, which is twice 
as strong as that obtained from /x — > eee decay. In particular, in the case 
of the Eq models [Q, these limits provide stringent constraints on the 
mass {Mz > 5TeV) and the mixing angle of Z and Z'. 

3.2.3. /5/?oj/ and {fi~ — e+) processes. We discuss now briefly the lepton 
number violating processes eqs. (H), (p!oD (see also ref. Il2|, |67|). 

• The oldest lepton violating process, intimately related to the nature of the 
neutrino, is the neutrinoless PPof decay 

{A, Z) ^ (A, Z ± 2) + + e^, + (A, Z) ^ {A, Z - 2) + e+ (106) 
which together with the allowed PP2u decay 



{A,Z) ^ (A,Z±2)+e^ + | |, e^ + {A,Z) (A, Z - 2) +e+ + z/, + 1/^ 

(107) 

are the only modes of some otherwise absolutely stable nuclei. One finds 
that |]12[, the life time of these processes is given by 



Ty2{oiy)=Ko./\n\'\ME\l (108) 



Tu2i2iy)=K2./\ME\i, (109) 
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where rj is the relevant lepton violating parameter and | ME \ the correspond- 
ing nuclear matrix elements. The quantities Kq^,, are functions of {A, Z) 



which can take the following values 



1.5 X lO^^y < Kq^ < 2.5 x lO^^y, (110) 



2.5 X lO^^y < < 3.3 X lO^^y (111) 

Thus, 7], which contains all information about the gauge models, should not 
be much smaller that 10~^ (the present experimental limit) in order l3(5oi, 
decay to be within the capabilities of present experiments. In the 5*0(10) 
model discussed above, the best value for r] is given in the case of light 



neutrinos r/^^ ~ 4 x 10~^. It is not difficult, however, to invent scenarios 
where (3j3ov decay is not far from the experimental capability of the near 
future experiments. 



In particular, there have been recently discussed cases [§8] with three light 
neutrinos predicting /3/3oy decay half-hfe in the range of the sensitivity of 
future experiments, induced by the exchange of majorana neutrinos with an 
effective neutrino mass ~ (0.1 — 1.0)eV . The interesting feature of these 
scenarios is that they also provide a solution to the solar neutrino problem 
and in some particular cases, they can also accommodate a solution to the 
atmospheric neutrino problem. 

• The other interesting lepton number violating process is the (/i~,e"'") con- 
version 

/x,- + (A^)^e+ + (A^-2)* (112) 
The half-life time is given by 

Ti/2(/x-,e+)=K^,/|nf|ME'|2^ (113) 

where, K'Q^{^l- , e+) = 2.2 x IQi^^yA^'^ / {Z'^^ff/ Z) or 5 x Wy <K'q^<6x lO^y, 
for nuclei ranging from ^^Ni to ^^C |^^. Even though this process is 10^° 
faster than its sister (e~, e"*"), unfortunatelly it must compete against ordinary 
muon capture /i^ — * z/^. Thus, one gets the branching ratio 
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^ r(/i-,z/^) ~ ^'^ ^ Z(1.62Z/A-0.62) 

where rj is the corresponding lepton violating parameter. Thus, for \vl\ ^ 
1^1 < 10~^, this process is unobservable. Even if transitions to all final nuclear 
states are considered [^], \ME'\^ ^ O.IZ^ = 100. The present experimental 
limit is R<9x 10"^^ j^ggg gg^^^ 2.3). 



3.2.4. Other lepton flavor violation mechanisms. We mention here some 
other mechanisms of lepton flavor violation which are also possible. Of par- 



ticular interest is the Bjorken- Weinberg ||6^ mechanism, based on a two loop 
contribution when additional Higgs bosons are present. The result is sensi- 
tive to the various unknown parameters (Higgs mass etc. [|7^), but it can in 
principle fall close to the experimental limit. 

Another possibility of lepton flavor violation arises, if we extend the scalar 
part of the theory with the introduction of a doubly charged singlet A 
already mentioned in the sect. 3.1.2. The most important effect of this 
scalar is in the — > 3e decay, which is mediated at the tree level [Q. One 
finds that the branching ratio is 

M.-3,4(£^)^(^)' (-) 
which, combined with the experimental limit, gives the following bound 

> 3.2 X lO^GeV^ (116) 



Before closing this section, we should recall that, lepton flavor violating pro- 
cesses play a crucial role in theories proposed to solve the hierarchy problem. 
One such example is given by the theories of dymamical symmetry breaking 
(Technicolor (TC) and extended Technicolor (ETC) [0.) Flavor chang- 
ing reactions, both in leptonic as well as in the quark sector, are found to 
be incompatible with the experimental limits in these theories. It has been 
shown recently |]74| that, fixed point or walking technicolor theories can solve 
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the problem of large flavor violation and bring it down to the experimentally 
allowed region. 

The most interesting theory, which can solve the hierarchy problem of cource, 
is the theory of Sypersymmetry. Lepton flavor violating reactions are always 
present in this case, receiving new contributions from the supersymmetric 
partners and renormalization effects. These, will be discussed in the next 
section. 



3.3. Supersymmetric extensions of the standard model 

One of the most important extensions of the standard theory is the minimal 
supersymmetric standard model (MSSM). Although there is no experimental 
evidence of supersymmetry as yet, it is a common belief that supersymmetry 
and supergravity play an important role in the theory of elementary particles. 
The main motivation for the incorporation of supersymmery in the funda- 
mental theory of interactions is the natural solution of the hierarchy problem 



I^. Moreover, supersymmetry seems to play important role in other issues 
of the unification program. SUSY models predict a longer life-time for pro- 
ton, they play crucial role in models with inflation, while they predict exact 
unification [|7S[ of the three gauge couplings consistent with the LEP-data. 
Furthermore, superstring theories [^, which appear today as the only can- 
didates for a consistent theory unifying all fundamental interactions, result 
to an effective superymmetric theory in low energy. 



Flavor changing neutral currents [^, |7^, |79|, ^ present one of the most 



important tests of all these low energy supersymmetric theories. In the fol- 
lowing subsections we are going to present a brief review of the effects of the 
new sources of flavor violation in the context of the above theories. 



3.3.1. Minimal Supersymmetric Standard Model. In the minimal su- 
persymmetric extension of the standard model one can write down the fol- 
lowing Yukawa couplings, 

W = XuQHU^ + XdQHD^ + XeLHE^' (117) 
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where Q, D , E , are the usual superfields which accommodate quarks and 
leptons. The potential is given by 



V = EI^I +m^/2El¥'d' + ^(W + W*) + S(|^|¥.. + c.c.) (118) 

where 7713/2 the gravitino mass and A, B are the scalar mass parameters 
depending on the details of the superymmetry breaking. Then, the s-lepton 
6 (S> 6 matrix in the basis (e, e*) takes the form 



m\i2l + 7nlme A*m\ 

A7ne ^1/2^ + 7n\7ne 



(119) 



where me is the usual lepton mass matrix and A = A + 2B. In this approx- 
imation, the lepton and s-lepton mixing mass matrices are similar which, 
in turn, implies that SlSg = 1. This means that, there is no lepton flavor 
violation induced. 

One can also show that, in this model there is no contribution from neu- 
tral intermediate particles, since at this level either they remain massless 
(neutrinos) or degenerate (s- neutrinos). 

The inclusion of the isosinglet right-handed neutrino, A^*" = {N^;Nc), can 
lead to additional terms in the superpotential of the form 

Wi = XnLHN^ + ^MnN^N^ (120) 

Now, the 6 ® 6 neutrino mass matrix can have both Dirac and isosinglet 
majorana mass terms. The corresponding 12 ® 12 neutral s-lepton mass 
matrix becomes analogous to that of the neutrinos. Even though 5*1 S*^ is 
non-diagonal, lepton flavor violating processes are suppressed due to the fact 



that s-neutrinos are degenerate |]78 



The above results are modified if one goes beyond the tree level and includes 
radiative corrections arising from the first term of the superpotential (eq. 
(|120D ) taking into account renormalization effects |^0|, ^ ^] and considering 
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Figure 3: Radiative contributions to the Li masses. 



radiative contributions to the scalar masses (Squarks, Sleptons etc,) at the 
one loop level. Indeed, let us assume that Li are the quark and lepton 
superfields and that the MSSM is extended with the inclusion of X, Y, Z 



additional superfields [pO]. Then, one gets the additional Yukawa couplings 



W = KjLiLjX + A'^L.rZ (121) 

The additional terms create the diagrams of fig. 3, which lead to radiative 
contributions of the type 

/\Ml (X A,,A,, + A;a; (122) 

where \ij is a 3 ® 3 matrix in generation space, while A^ is a column vector. 
Now, if the fields X, F, Z are light compared to the Plank scale, then the 
above corrections are proportional to a large logarithmic factor ln{Mpi/ Mx,y.z) 
and the contributions to the scalar masses are significant. The origin of fiavor 
violation here lies in the fact that, due to these contributions, it is no longer 
possible to diagonalize fermion and s-fermion mass matrices simultaneously. 

Let us consider the above effects in the case of s-leptons. In the presence 
of the right handed neutrino, the mass matrices discussed in the previous 
section, receive contributions which modify the tree level results as follows 

mgrnt = m\i2l + rnem\ + cmjj^m\^ (123) 

where rrty^ is the Dirac type neutrino mass matrix as before. The coefficient 
c = c{t) is scale dependent {t = Infj,) and includes the running from the 
Plank scale down to the scale where the right-handed neutrino acquires a 
mass M]\r ~ 0{Mgut), thus c is proportional to ln{Mpi/ Mgut)- Because 
the dependence of the corrections on the scale is logarithmic, and the range 
Mpi — Mgut is not very large, the corrections inducing the flavor mixing in 
the charged s-lepton mass matrix can be written as 
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(Am; 



3m2/2 + A 



rln- 



PI 



_{27Tvsinf3)^"~ Mgut 

where v = 2A6GeV and P = tan-^{< H > / < H > 
ized Dirac neutrino mass matrix at the GUT scale. 



(124) 



m' 



is the diagonal- 



The most natural candidate models, that the above analysis can apply, are 

|, Some particular 



the models derived from the superstring [p3 
cases will be considered in the next section. 



3.3.2. Flavor violation in Superstring Models. One of the most promis- 



ing avenues beyond local supersymmetry is the theory of superstrings 
Indeed, superstring theory is the best candidate for a consistent unification 
of all fundamental forces including gravity. On the other hand. Grand Uni- 
fied Theories are incorporated naturally in the superstring scenario. Realistic 
string models suggest that, all gauge interactions should unify within a simple 
gauge group, with a common gauge coupling gstring at the String unification 



scale M string, which is found relatively high and close to the Plank scale p7 



Mstrrng ~ gstr^ng X Mpi 5. X lO^^GeV (125) 

On the other hand, renormalization group calculations have indicated that 
minimal supersymmetric Grand Unified Theories (SUSY-GUTs) are in agree- 
ment with the precision LEP data when the SUSY-GUT scale Mq is taken 
close to Mg ~ lO^^GeV [^. Thus, the minimal supersymmetric standard 
model (MSSM) cannot probably be derived directly at the string scale; the 
above descrepancy between the two scales, would rather suggest that, either 
the MSSM should be obtained through the spontaneous breaking of some 
intermediate GUT- like gauge group which breaks at the scale Mq, or some 
extra matter fields are needed to modify the gauge coupling running. 

Thus, string motivated non-minimal extensions of the MSSM, contain ad- 
ditional Yukawa interactions, which may in principle lead to interesting en- 
hancement of flavor changing neutral processes. In particular, the radiatively 
induced lepton flavor violations, discussed in the previous section, occur nat- 
urally in this kind of models. 
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In order to be specific and give some quantitative results, we will concentrate 
on some realistic string constructions discussed extensively in the literature. 
These attempts have been made in the context of the free fermionic for- 
mulation of the four dimensional superstrings, and led to the construction 
of three types of models. Two of them are characterized by an intermediate 
GUT scale based on the symmetries SU{5)®U{1) Q and SU{4)^0{4) [H], 
while there is the third type of models where the original string sym- 
metry breaks down to the standard model symmetry times some additional 
U{1) factors. 



In the case of flipped SU(5) model the superpotential of the minimal 



model (assuming invariance under H —>■ —H) reads 



W = X^FiFjh + X'iFJjh + X'ifil^^h + X^HHh 

+ X^HHh + Xl^F.H^a + Xjhh^o + fiab^a^b (126) 

The Fi + fi + II (i = 1,2,3) are the three generations of 10, 5 and singlet 
representations of SU(5) that accommodate the light matter particles of the 
Standard Model, H and H are 10 and 10 Higgs representations, h and h 
are 5 and 5 Higgs representations, and the 4>o,(f)a (a = 1,2,3) are auxiliary 
singlet fields. The first 3 terms in the superpotential eq. ( |126|) give masses 
to the charge 2/3 quarks Ui, charge -1/3 quarks rfj and charged leptons Zj, 
respectively. The next two terms split the light Higgs doublets from their 
heavy colour triplet partners in a natural way. The sixth term provides a 
large element in the see-saw neutrino mass matrix, and the term XjhfKpo 
gives the traditional Higgs mixing parameter. Under the same assumptions 
discussed in the previous sections, the induced mixing in the s-lepton mass 
matrix, due to renormalization group running in the range Mpi — Mqut, is 



B|,|8|] 



{Smj) 



where, due to the fact that up-quark and Dirac neutrino masses arise from 
the same superpotential term, we have substituted m^^ with m^, which is 
the diagonal up-quark mass matrix at the GUT scale. 



37 



The diagonalizing unitary matrix V = SlS^, is unknown, due to the fact that 
nie and mass matrices are unrelated in the flipped SU{5). Thus, it is 
possible that, the mixing may enhance the flavor changing reactions in this 
model. 

In the case of the SU{A) ® 0(4) = SU{A) O SU{2)l ® SU{2)ji string model, 
the lepton and KM-mixing matrices are related due to the fact that leptons 
and down quarks receive masses from the same superpotential term. The 
superpotential of the model under a Z2 symmetry H ^ —H can be written 
as follows 



W = XiFLFRh + X2FRH<Pi + X^HHD + XiHHD 

+XAFrD + X^FlFlD + X^(t)^hh + X^(t)^ ( 128) 

The phenomenological implications of the above superpotential terms have 
been discussed extensively in previous works [0 . 



The minimal supersymmetric version of the model includes three generations 
of quarks and leptons, which are accommodated m. Fl + Fr = (4, 2, 1) + 
(4,1,2) representations of the SU{4:) SU{2)l ® SU{2)r symmetry. One 
needs at least one pair of H + H = (4,1,2) + (4,1,2) Higgses to realize 
the first symmetry breaking SU{4:) x SU{2)r — > SU{?>) x U{1)b-l, and one 
h = (1, 2, 2) higgs field to provide the two Weinberg-Salam doublets after the 
first symmetry breaking. In addition, a sextet field D = (6, 1, 1) is needed to 
produce a pair of colored triplets 3 + 3, which are going to combine with the 
uneaten dj^ and d'jj, to form two superheavy massive states and avoid fast 
proton decay. 

Thus, in the case of the minimal version of this model, one obtains the mass 
relations rriu = m^^ and me = rud at the GUT scale, which in turn imply 
that the lepton and KM mixing matrices are also equal at Mqut- This fact 
leads to a definite prediction for the above processes. Due to the small mixing 
angles, however, there is a significant suppresion of flavor violating processes. 



One novel feature of both models is the generalized see-saw mechanism [84 



85, 91 1 . Indeed, three of the singlet fields introduced, are used to realize the 
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see-saw mechanism together with the left and right handed neutrinos. The 
see-saw matrix in the basis (z/j, Nf, 0^) takes the form 



/ 




(129) 



V M^c^^ 

where it is understood that, all entries in eq. (|129|) represent 3x3 matrices. 



The m^^, M^c^^, submatrices arise from the trilinear superpotential terms 
of the above models. The contribution M^"*^ is a majorana type mass matrix 
for the right handed neutrino and usually arises from some higher order non- 
renormalizable terms. In principle, there are many arbitrary parameters in 



eq. ( |129| ), but within the context of some recently proposed fermion mass 
matrix Ansatze, as well as under some natural assumptions |^2|, |93[, it is 
possible to reduce the arbitrariness and obtain definite predictions. In what 
follows, we are going to use the results of a fermion mass matrix Ansatz 
proposed in the context of the above string derived models to estimate the 
renormalization effects on the various flavor violating processes. 



3.3.3. Flavor violating processes, i) We start again with the simplest 
decay 67. The diagrams are presented in fig. 4. We can parametrize 
the amplitude with the functions /mij/ei in the same way as in sect. 3.1. 
Here they are given by 

1 77i^ 

fMl = -fEl = --m^^f{x) (130) 

2 

where fj is the corresponding flavor violating quantity defined as follows 



V = - '^^"J':"' (131) 



and {Srhgg*) 12 is determined by the analysis of the previous section while rha 
is the mass of the haviest sparticle circulating in the loop. If = m| for 
example, it is given by p4|, BSl 
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Figure 4: jj, ^ e'y in supersymmetric theories. 



ml = ml/^ + Ceit)ml/^ (132) 

where Cgit) is a function of the scale t = Infi, and for fi ~ mw, Cg ~ 0.50. 
The function /(x) of eq. ( |130| ), depends on the ratio x = m^/rha, being 
the photino mass, and is given by 



1 TTi? 

f(x) = — -{1 + 2x^ + _ 6x - QxHnx}, x = —!r (133) 

12(1 — x)* 

Now, the branching ratio -Re7 takes the form 

i?e, - - (^^^^2)2 - 1^1 ^0 (134) 



where Rq is defined through 



As an apphcation, we use the results of ref. [^] for the leptonic mixing angles. 
Taking as an example the initial condition at the GUT scale, m^i2 = m\i2i 
we can obtain the bound m3/2 > 1'bGeV . 

a) The decay yU — ^> 3e is treated similarly. The corresponding diagrams are 
shown in fig. 5. We find 



(136) 

with 
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Figure 5: Diagrams for the /x ^ 3e process in supersymmetric theories. 



g(x) = — --{2 - Ux^ + ISx^ -9x + Qx'^lnx} (137) 
36(1 — xp 

1 771- 

fbi^) = - + '^^ + + '2)lnx}, X = ^ (138) 

8(1 — x)^ 

Comparing — > 67 with /i ^ 3e we get 

— ^ ~ — {16 } - 6.4 X 10"^ 139 

Re^ 247r*- me 3 6 -' ^ ^ 

in) Let us now discuss the {lJ> — e) conversion. The amphtude for this process 
is given by 

M^I^^'^C} (140) 

where the first term corresponds to the photonic and the second to the non- 
photonic contributions arising from the box diagrams (see fig. 6) and C = 
"^3/2/"^l- We find that 

= m(pi)(/mi + 75fEi)io-^''— + (/eo + 75/Mo) Y (gxiy - 

(141) 

= jV7A^^^iV, N = Nucleon (142) 
For the box diagrams we obtain 

J('2) = «(Pi)7'^(/v + /475)«(Pm) (143) 
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Figure 6: Diagrams for the {/i — e ) conversion in supersymmetric theo- 
ries. 



Jx^ = Nlxl [(3 + Pfvrs) - ifv + fAf3r^)-f,] N (144) 



where (3 — (3q/(3i, with 



Furthermore, 



9 9 mi 

a 



1 771^ 

/eo = -/mo = --^m^gi^)^- (147) 



3/2 

/y = -/A = -^f,a'h{x)-4- (148) 



2 

"^3/2 



It is quite hard to write down general expressions containing both photonic 
as well as non-photonic diagrams (sec sect. 4.1). For the coherent process, 
however, one can write down the (/^ — e) conversion rate as follows 

1 1 ~ |2 1 ~ |2 1 

~ {G m?Y ~^ ~^ 2'^^^ ' ~^ ' ~^ "^^^ ~^ 2 '^■^'^ 

(149) 

where 

^-('l + ^^^V (150) 
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and 



Fz,N{.(f') are the nuclear form factors to be discussed later and /(A, Z) is the 
Primakoff's function (see eq. (p9|)). Under some plausible approximations, 
we can write the (/i — e) branching ratio as 

Assuming that the photino is the lightest supersymmetric particle, we can 
take a; ^ 1, to obtain 

111 

/ = — , g = —, fb = - (153) 

12 18 8 ^ ^ 

Thus, the two terms from the photonic contribution tend to cancel and the 
box diagram dominates. Comparing the (/i~, e~) conversion with the — 67 
decay, we get 

ReN ~ ^ (^3 + Xj ^P'^- ^^^^ 

Analytic results of the abovebranching ratios are presented in sect. 5.3. 



4. EXPRESSIONS FOR THE BRANCHING RA- 
TIO OF (/i-,e-) 

In this section we will discuss the construction at nuclear level of the effective 
operators which are responsible for the semileptonic process [i~ (A, Z) — »• 
e~(A, Z)* and derive the expressions for the branching ratios of this process. 
If we assume that the A nucleons of the nucleus (A,Z) interact individually 
with the muon field (impulse approximation), then the needed nuclear matrix 
elements can be obtained from an effective Hamiltonian f2 which arises from 
that of the free particles, as it is described in the next section. 



4.1. The effective Hamiltonian of (/i ,e ) conversion 
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As we have seen in sect. 2, the effective amphtude for (/x , e ) conversion 
is given by eq. ( p.4CI| ). It is not easy to separate the dependence on the 



nuclear physics from the leptonic form factors. This can only be done for 
the coherent mode. For the general case however, we will discuss separately 
the photonic and the non-photonic contributions. For compactness of our 
notation, we will write the hadronic currents given in eqs. (|142|) , ( 144 ) in 
the general form 



where 



J^? = N^x [^y(3 + Pfvrs) - Ufv + fAPr^h^] N (155) 



gv = 7:, Qa = 0, fv = 1, P = 3 {photonic case) (156) 
6 



dv = dA = 2' fv = 1) fA = 1-24 {non — photonic case) (157) 

The parameters /3, Jeo, /ei, Imo, Imi, Ia, fv and C depend, of cource, on the 
model. Thus, in the case of the neutrino mediated processes in left handed 
theories one finds 



_ Ci^m^ /3_/30l ^_/25 {light neutrinos) 1 ^ip-o\ 
^ ~ ^2~' ^° ~ \ 1 J ' '^^ ~ I 5/6 {heavy neutrinos) j ^ ^ 



i.e. /? = 5/6. Furthermore 

fv = Pofi, fA = M2 (159) 

The quantities /i, /2, depend on the specific gauge model and are given 

in the literature |T2|. For the supersymmetric models, the corresponding 



expressions have already been given in sect. 3.2.2. 

At the nuclear level the relativistic expression of the hadronic current is not 
needed. Using the standard non-relativistic limits the relevant nuclear matrix 
elements involve the operators 
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where the summation is over all nucleons and Pe is the momentum of the 
outgoing lepton. The factor l/-\/3 is introduced to make the two matrix 
elements in the total rate expression equal and compensated by suitable 
factors elsewhere in the expression. In eq. ( |160D 



e = fvlfA = 1/1.24 (161) 

Assuming that the kinetic energy of the final nucleus is negligible and taking 
me ~ 0, we can write the magnitude of the momentum pe of the outgoing 
electron approximately as 



Pe 



m,. 



Cfo — {Ef — E, 



gsj 



(162) 



where Ef, Egg are the energies of the final and ground state of the nucleus, 
respectively, is the muon mass and eb is the binding energy. For coherent 
processes Ef = Ei and since is relatively small (the biggest value occurs 

0.53/m~^. For 



in lead region where ej, ~ lOMeV) we have p, 
incoherent processes Ef ^ Ei and for a sum-rule approach one can assume 
a mean energy |]63| , |6^ for the emitted electron corresponding to a mean 



excitation energy of the nucleus E. In this case p, 



E-eb. 



By expanding the exponential of the operators in eq. (|160|) in terms of spher- 
ical Bessel functions ji{x), we obtain the multipole expansion of the (/i~, e~) 
operator, i.e. the following two types of operators T*^'-'^)"' 



i=l 



for the vector part, and 



(163) 



T, 



M 



9a\Iy II(^ + /^^3i)jz(?^i 



M 



(164) 
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for the axial vector part. In eqs. ( |163| ), ( |164| ) q represents the magnitude of 
the momentum transferee! q to the nucleus during the e~) process. In a 
good appromimation g ~ Pe and thus, q is given by the energy conservation 
eq. (^). 



4.2. Expressions for the branching ratio of (/i , e ) conversion 

The probability density for converting the bound /i~ of a muonic atom to an 
e~ with momentum pe is given by the Fermi's golden rule 



r„; = I/*. (I;) \<S\a\i,,>\'-^ (i65) 

where | i, /i > is the initial state of the system: nucleus (A, Z) + and 
I / > the final state of the system: nucleus (A, Z)* + e~ . The factor {pe/fn^Y 
involves the density of the final states appropriate for normalization of the 
wave packet representing the outgoing e~ and Pe is the unit vector in the 
direction of the electron momentum. The quantity S depends on the gauge 
model (in the photonic case, for example, it coinsides with of eq. (^7])). 
The effect of nuclear recoil has been neglected. In the above expression 



\< f \ Vt\i,n >P=|< / I I >P +^\<f \ \ i,^^>\'^ (166) 

Of special interest, as we shall extensively discuss below, is the partial rate 
for the gs gs transitions i.e. the coherent rate. It has been previously 
assumed that, the Is muon wave function varies very little inside the light 
and medium nuclei [p6|, i.e. the following approximation has been used 



where 



<f\n\ >|^=< $1, >"|< f\n\i>\' (167) 



2_ Id^x\<l>,{^)\^p{^) 

<'^'^>= — JdF^) — ^^^^^ 

In the latter definition, $^(x) is the muon wave function and p(x) is the 
nuclear charge density. To a good approximation it has been found 
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< $1, >2= ^ (169) 

TT Z 

{a is the fine structure constant) i.e. the deviation from the behaviour of 
the wave function at the origin has been taken into account by Zejf. The 



above approximation was first used for the (/i , v^) ^ and afterwards 



in analogy to muon capture in the (/x , e ) process 0. Recently 



103|| , it had been found that eq. (|167|) is not very accurate in ordinary 



muon capture. We will therefore present results, both, with and without this 
approximation. 

The branching ratio of the total (/x^, e^) conversion rate divided by the total 
muon capture rate, assuming the approximation of eq. (|167|) for both (/x, e) 
and /i, v) processes and only photonic or non-photonic mechanisms, takes the 
simple form 

ReN = = PI (170) 

r(/i, z/) 

where the quantity p contains all the nuclear dependence of ReN and 7 con- 
tains the elementary particle physics. We mention that ReN is the quantity 
provided by experiments. Obviously, the effect of the approximation eq. 
( |167|) on the branching ratio is expected to be negligible. 

Another interesting quantity in the study of the /i — e process is the ratio of 
the coherent (/i~,e~) rate, Veoh = ^gs-,gs, to the total (/i~,e~) rate, Ttot = 
J2f^i^f for all final states | / >, i.e. 

V = (171) 
This will be discussed in detail below. 



4.2.1. Coherent (/x , e ) conversion. 

In the case of the coherent process, i.e. ground state to ground state (0"*" — > 
O"*") transitions only the vector component of eq. ( |166| ) contributes and one 
obtains 
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<f\no\t,fi> = gv (3 + /y/3)F(g2) 



(172) 



where F{q^) is the matrix element involving the ground state 



d'x{p,{^) + ^_Mp„(x)}e-''-<|.^(x) 



(173) 



3 + fvfi' 

(pp(x), Pn{^) represent the proton, neutron densities normalized to Z and N, 
respectively). F{q^) can also be written as 



-J + JvP 



(174) 



where 



F.J 



d^x pp,n{x) e <l>^(x) 



(175) 



In eqs. (|174|) , (|175|) we have not made use of the eq. ( |167| ). Thus, the coherent 
rate can be obtained from eq. (|172|) by first calculating -F(g^) from eq. ( |173| ) 
for a given muon wave function and a given nuclear density distribution. 
In ref. the muon wave function was obtained by solving numerically 



the Schrodinger equation taking into consideration the effects of vacuum 
polarization and finite nuclear size. Weinberg and Feinberg [0 used for the 
coherent rate the expression of eq. ( |167|) and estimated the quantity Z^ff 
as in ref. 
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Shanker using Fermi nuclear distribution included an 
additional correction interference term. Using eqs. ( |1671 ) and ( |169| ), for the 
coherent rate one obtains 



\F{q 



2m2 



IT 



'eff 



ZFz{q^ 



3-/y/5 



2m2 



with Fz [Fn) the proton (neutron) nuclear form factors 



NFNiq')\ 



(176) 



1 



1 



c/=^xp„(x)e-''i-" (177) 



These nuclear form factors can be calculated by using various models i.e. 
shell model ||100|| , quasi-particle RPA ||101|| etc., or can be obtained from 
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experimental data whenever possible ||102|| . The branching ratio in this ap- 
proximation takes the form 



eN 



7.2 



2^2 



9v 



1 + 



3 + /y/? 



Ivh. 



(178) 



with 7p/i defined in eq. ( |151|) . The validity of the approximation of eq. ( |167| ) 
will be discussed in sect. 5.3. The parameters F/v, -Fz, Fp, Fn and Z^ff 
for various nuclear systems appear in table 1. The corresponding widths in 
arbitrary units are given in table 2. 



Table 1. Parameters needed for calculations of the coherent (yU~,e^) conver- 
sion matrix elements with: i) the exact muon wavefunction, Fp and Fn of eq. 
( \17d{) , and ii) two-parameter Fermi distribution, proton and neutron form 
factors, Fz and F^- The electron momentum pe and the effective charge 
Zeff of a set of nuclei covering the hole periodic table are also shown. 



A 


Z 




Fnlfm-^/'] 


Fz 


Fn 


Pe{MeV/c) 


^eff 


12 


6 


0.86 10^2 


0.86 10-2 


0.77 


0.77 


105.067 


5.74 


24 


12 


0.37 10-1 


0.36 10-1 


0.65 


0.64 


105.017 


10.81 


27 


13 


0.45 10-1 


0.45 10-1 


0.66 


0.62 


104.976 


11.62 


32 


16 


0.66 10-1 


0.63 10-1 


0.62 


0.59 


104.781 


13.81 


40 


20 


0.99 10-1 


0.92 10-1 


0.58 


0.54 


104.449 


16.47 


44 


20 


0.97 10-1 


0.11 10° 


0.57 


0.55 


104.464 


16.43 


48 


22 


0.11 10° 


0.13 10° 


0.55 


0.52 


104.258 


17.61 


63 


29 


0.17 10° 


0.19 10° 


0.49 


0.46 


103.474 


21.22 


90 


40 


0.26 10° 


0.29 10° 


0.42 


0.38 


101.951 


25.69 


112 


48 


0.29 10° 


0.36 10° 


0.36 


0.33 


100.749 


27.86 


208 


82 


0.42 10° 


0.57 10° 


0.25 


0.22 


95.125 


33.81 


238 


92 


0.40 10° 


0.57 10° 


0.20 


0.18 


93.591 


34.36 



Table 2. Coherent widths (in arbitrary units) for the photonic and non pho- 
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tonic mechanisms: (i) with the exact muon wave function, $^(x) (ii) with 
< $j^(x) > in the approximation of Z^ff. The ratios of the width in case (i) 
to the width of ordinary muon capture with arbitrary normalization are also 
given. 



non-photonic mechanism {(3 — 5/6 ) photonic mechanism (/3 = 3) 







With 


With 




With 


With 










A 


Z 




< $.(x) 


> 


Ratio 






< ^/.(x) > 


Ratio 


12 


6 


0.52 10-^ 


0.51 10- 


-4 


0.14 


0.21 10" 


-4 


0.21 10- 


-4 


0.058 


24 


12 


0.91 10-3 


0.90 10- 


-3 


0.19 


0.38 10- 


-3 


0.37 10- 


-3 


0.078 


27 


13 


0.14 10-2 


0.14 10- 


^2 


0.20 


0.56 10- 


-3 


0.55 10- 


-3 


0.083 


32 


16 


0.29 10-2 


0.28 10- 


-2 


0.22 


0.12 10- 


-2 


0.12 10- 


-2 


0.093 


40 


20 


0.64 10-2 


0.62 10- 


-2 


0.26 


0.27 10- 


-2 


0.26 10- 


-2 


0.110 


44 


20 


0.72 10-2 


0.69 10- 


-2 


0.40 


0.26 10- 


-2 


0.25 10" 


-2 


0.150 


48 


22 


0.94 10-2 


0.91 10- 


-2 


0.36 


0.35 10- 


-2 


0.34 10- 


-2 


0.140 


63 


29 


0.21 10-^ 


0.19 10- 


-1 


0.36 


0.78 10- 


-2 


0.72 10- 


-2 


0.140 


90 


40 


0.47 10-1 


0.42 10- 


-1 


0.54 


0.17 10- 


-1 


0.16 10- 


-1 


0.200 


112 


48 


0.62 10-1 


0.54 10- 


-1 


0.62 


0.22 10- 


-1 


0.19 10- 


-1 


0.220 


208 


82 


0.13 10° 


0.89 10- 


-1 


0.98 


0.41 10- 


-1 


0.29 10- 


-1 


0.310 


238 


92 


0.12 10° 


0.73 10- 


-1 


0.92 


0.35 10" 


-1 


0.22 10- 


-1 


0.280 



4.2.2. Incoherent {fj,~,e~) conversion. As we have mentioned in sect. 
2.3, from the experimental point of view, the coherent contribution to the 

{fi~, e~) conversion is the most interesting. It is, however, important to know 
what fraction of the total rate goes into the coherent mode and how this varies 
with A and Z. In this section we will evaluate the incoherent contribution 
and will compare it to that of the coherent channel. 

For the calculation of the contributions leading to the exited states many 
methods exist which involve various approximations. Those which are based 
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on the aproximation inserted by eq. ( |167| ) will be discussed in detail in sect. 
5. In this section we will elaborate on a method recently developed by Chiang 
et ai, and which uses the exact muon wave function in eq. ( |165| ). This 



method, known as "nuclear matter mapped into nuclei by means of a local 
density approximation" , is accurate and leads to a good reproduction of the 
/i~ capture rates over the periodic table once the proper renormalization of 
the weak currents is considered. 

In the ordinary muon capture only the incoherent channel is open. In this 
process the total rate T^c is given by | p.03| 



r^e = /rfV|<i>^(x)|2 r(p„(x),pp(x)) (179) 

where r(p„, pp) is the width in an infinite slab of nuclear matter with neutron 
and proton densities Pn,Pp, respectively. This width can be obtained by 
means of the Lindhard function as 

r(p„,pp) = -2| ^n,^SE|r|2 im%,nip, - pu) (iso) 

where T is the transition matrix, i.e. the amplitude associated with the 
elementary process fi^p nv^ and mi,Ei are the masses, energies of the 
particles involved in this reaction. The Lindhard function U corresponds to 
ph excitations of the p-n type and is given by (see e.g. ref. ||104| , |106|| for its 
usage) 

fl (n\-0 ( "^^P ni(p)[l-n2(q + p)] 

with 7^1 (q') and n2(q) the integral (0 or 1) occupation probabihties of the 
neutron and proton, respectively. 

In an analogous formalism, the incoherent (p~, e~) conversion rate, rj„(;(p — > 
e~) = J2f^gs Tj-*/, can be expressed as follows 

r_(p-A-.e-x) = -2 /rf3x|<|.,(x)p/|%n,^^^ 



J (27r)3 

X [SS|Tp(p"]9 e'p)ImUp^p{p^ - Pe) 

+ T.J:\T\'^{p~n e'n)ImUn,n{Pfi - Pe)] (182) 
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The two terms in the brackets of eq. ( |182| ) result from the character of the 



(/i , e ) operator to be of charge conserving type. 

One can separate the nuclear dependence from the dependence on the el- 
ementary sector in the two processes, (/i~,z/^) and (/i~,e~), by factorizing 
outside the integrals in eqs. ( |179D , ( |182D an average value of the quantity 



SS|Tp. For the ratio of the incoherent (/i^,e~) rate divided by the total 
muon capture rate, since SS|Tp for the first process is proportional to -E'^-E'e 
while for the ordinary muon capture it is proportional to E^Ey, we get 

^ ^ ^inc ^ mpmJ^{E^Ee)~^T^\T\^{^-p e-p)]a^ ^ 
r^c mnmy[{E^Ey)'-'^J:j:\T\^{fi-p ^ nu^)]av ^ 

The quantities Gp, Gn are smooth functions of the momenta and contain the 
nuclear dependence of the ratio R. They are defined as 



IJ. Pe) r'\QA\ 









/ d^peImUp^p{p^ - Pe) 






P(2vr)- 


"3 / d^pJmUp^n{p^i - Pu) 



r (7 m- Z Jd^x\<!>f,{^)W27TfJd^pJmUn,n{P,^-Pe) . . 

' ^ NJd^x\^^{^m2nrjd^pJmUp,^ip,-p,) ^ ^ 

Note that R differs from the branching ratio ReN, because R does not include 
the coherent part of the (/i~,e~) process. The calculation of R in various 
nuclei is based upon Gp and Gn, since the |Tp at the elementary level is in 
most models independent of the nuclear parameters (see ref. ||103|| for values 
of |Tp in ordinary muon capture). As an example, we give here the ratio R 



in the non-photonic case |T^ involving the box diagrams of sect. 3.3.3 which 
is 

All the nuclear information is contained in the function G{A, Z) defined as 
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G{A,Z) = Gp{Z,N) 



Pv + 3/1 



+ GniZ,N) 



N 




(187) 



We should mention that the quantities Gp, Gn contain Pauh blocking cor- 
rections and for the ordinary /i~ capture reaction the Q value, which is sig- 
nificant for light nuclei, is also considered in the argument of the Lindhard 



function in eq. ( |181| ). In heavy nuclei the Q value is approximately given in 
terms of the difference between the neutron and proton Fermi energies and 
no other correction is needed. 

As a summary we would say that, the use of an exact muon wave function 
for the calculation of the coherent rate, gives more accurate results in the 
medium and heavy nuclei region. It involves an approximation (local density 
approximation) for the sum of the non-coherent contributions. It uses neither 
closure nor explicit summation over all final states. It makes a summation 
over the continuum of excited states in a local Fermi sea. Its accuracy is 
tied to the number of states one can excite; the larger the better. The muon 
mass provides the energy for the excitation of such states. Furthermore, the 
method is quite simple to apply in actual calculations. 



5. EVALUATION OF THE NUCLEAR MATRIX ELEMENTS 

In this section we will discuss the methods of calculating the nuclear matrix 
elements needed for the partial and total (/i~,e~) conversion rates. It is 
known that the main feature of the (photonic) /i-e conversion rates is the 
strong Z-dependence (see below eq. ( |188| ) ) accounted for by assuming that all 
protons can interact independently with the muon and that this interaction 
is proportional to the muon density in the Is atomic orbit at the position of 
the nucleus. 

In this case, as we have explained in sect. 4, we can assume that, the prob- 
ability density for converting a bound muon in the Is orbit to an electron 
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of momentum pe, is approximately analogous to the muon average proba- 
bility density over the nucleus, < $is and the nuclear matrix elements 
Mj^ between the initial and final nuclear states. Thus, in order to find the 
partial rate Ti^f or the total rate F = J2f^i-*fj one needs to evaluate the 
contribution of the nuclear matrix elements | < f\^l\i > P for the states |/ >. 



5.1. The coherent /i — e conversion matrix elements 

In the approximation in wich one factorizes outside the integral of eq. ( |173| ) 
the average value of the muon wave function (see eq. ( |167|) ), the 

nuclear dependence of the rate in the coherent process, analogous to the 
matrix element M^^^g^ , is written in terms of the proton and neutron elastic 
nuclear form factors, Fz and F/v (see eq. ( |176|) , respectively, as 



M. 



gs^gs 



9H3 + fv(3rZ''FUq' 



fv(3 N F^jq^) 
fvP Z Fziq") 



(188) 



We note that in the photonic case (/y = l,/9 = 3) the neutron contribution 
vanishes and the nuclear matrix element becomes Z'^Fz{q^Y ■ We also men- 
tion that the nuclear form factors for the coherent — e process are calculated 



at q 



m, 



.534/m- 



The nuclear form factors, Fz{q^) and F^i^q^), have been calculated by using 
various models i.e. Fermi distribution |9^, ^ shell model ||100|| , quasi-particle 
RPA ||101|| . In the framework of shell model for spherical nuclei it was found 
that Fz{q^) and F]^{q^) can be cast in tractable analytical forms containing 



fractional occupation probabilities ||107|| , which take into account a significant 
part of the nucleon-nucleon correlations. Recently [p.06|| , the neutron form 
factors Fis[{q^) have been also extracted from the analysis of pionic atoms 
by means of a two-parameter Fermi distribution for heavy and very heavy 
nuclei and a harmonic oscillator density for light nuclei. The form factor Fz 



can be obtained from the available electron scattering data [102 . 



In the context of quasi-particle RPA one can calculate the coherent /i — e 
matrix elements of eq. ( |188| ), by using as ground state an uncorrelated or 
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correlated vacuum. In the first case the nuclear form factors, Fz and F/v 
take the form 



Fz{q^) = ^E(2J + 1) <j||jo(gr)||j> (V/)' (189) 



Mq') = ^E(2J + 1) <j\\jo{qr)\\j> {vff (190) 



The quantities Vf , Vf^ are the amplitudes for the proton, neutron single 
particle states to be occupied. Their values are determined by solving the 
known BCS equations iteratively ||108||, lie between one and zero and differ 



from those involved in the independent particle shell model (0 or 1). This 
is due to the consideration of pairing correlations in the RPA ground state, 
which deforms the Fermi surface of the nucleus, a picture described by the 
fractional occupation probabilities. The corresponding shell model form fac- 



tors Fz with fractional occupation probabilities, have been determined ||107 
by fitting to the electron scattering data. 

The accurate determination of the RPA ground state is of great importance 
for the exact calculation of the coherent and the total (/i~,e~) rate. The 
g.s. wave function provides the gs gs transitions, which are the dominant 
channel of the (/i~, e~) process and the total rate in the sum rule approach. 
A usual correction inserted in the ground state is the g.s. correlations | |109| , 



110| , |111| , |112| , |113|| , which can be included in the ground state by defining the 



correlated QRPA vacuum | > in terms of the uncorrelated vacuum | >. 
By using the Thouless theorem the correlated vacuum | > can be written 
as 



> = A^oe^^l > (191) 



where the operator 



S^ = IY: ,^C^l^^^A^a,JM)A^{b,JM) (192) 

'^ab,T,.JM ^'J -r ^ 

The operators A~^{a,JM) denote the two quasi-particle (or pair) creation 
operators, in the angular momentum coupled representation. The indices a 
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and b, denote proton (r = 1) or neutron (r = —1) configurations coupled to 
J, i.e. a = {j2 > ji) (and similarly for b), with jj running over the single 
particle states of the chosen model space: ji = {ni,li,ji). The correlation 
matrix C (symmetric matrix) is constructed for each multipole field A from 
the X and Y matrices i.e. from the RPA amplitudes for forward and backward 



excitation. A usual approximation for C is the following jlll 



a 



(A) 
ab 



(fW[xW]"')^^ (193) 

In eq. ( p.91|) , Nq is the normalization constant, which by keeping terms of 
first order in the correlation matrix C is given by 



1 



2 / ' ab ab 

ab,\,T 



V{A,r) q{X,t) 



(194) 



By using as ground state the correlated RPA vacuum of eq. (|191D , the 
coherent rate matrix elements take the form 



<o|r|o>=A/'o2<o|r|o> (195) 

which means that the correlated matrix elements are a rescaling of the uncor- 
related ones (see a similar expression in sect. 5.2.2 for the total rate matrix 
elements) . 



5.2. Total (/i , e ) Conversion branching ratios 

To find the total (/i^,e^) conversion rate, one need evaluate the matrix 
elements for both the vector and axial vector operators (see eq. ( p.66| ) ) and 
for all the final nuclear states | / > i.e. the quantities 

(q/ is the unit vector in the direction of the momentum transfer q^). Then 
the total (yU~, e~) rate matrix elements are given by 

Mtot = Sv + SSa (197) 
For the calculation of 5*^ and 5*^, one can use the following general methods: 
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1) Summing over partial rates: 



With this method we construct exphcitly the final nuclear states | / > in the 
context of a nuclear model e.g. shell model, random phase approximation, 
etc. The total e~) rate can be obtained by summing the partial rates for 
all possible excited states in a chosen model space. By using the multipole 
expansion of the (/i~,e~) operators eqs. (|163|) and (|164|) , the total rate 



matrix elements M^^^ can be written as 



tot 



^ fexc 



( las 



l,J 



II f II >| + 
f(''")-^ II >| 



198) 



LJ 



(s=0 for the vector operator and s=l for the axial vector one). The first term 
in the brackets of eq. ( |198| ) contains the contribution coming from all the 
excited states | fexc > of the final nucleus (incoherent rate) and the second 
term contains the contributions coming from the ground state to ground state 
channel (coherent rate). 

2) Closure approximation: 

It is well known that for the description of the total strengths in many pro- 
cesses the sum rule techniques are very useful |^, ^ |114| , p.l5|| . In such 



^gs 

This 



an approach the contribution of each final state | / > to the total rate is ap- 
proximately taken into account without constructing this final state explic- 
itly. One assumes a mean excitation energy of the nucleus E =< Ef > —E, 
and uses closure for the final states | / > i.e. J2f\f >< /I = 1- 
approximation requires the evaluation of one and two-body matrix elements 
involving only the initial (ground) state. This way the explicit calculation of 
the final states | / > is avoided. 

The mean excitation energy E of the nucleus, involved in the needed matrix 
elements is defined as I1B3 



E 



M 



2 , , - , , ,2 

\<f\TJ\gs>\ 



(199) 
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The numerator of this definition represents the energy weighted sum rule and 
the denomenator the non-energy weighted sum rule. Though E is defined 
in analogy with the ordinary muon capture reaction, the value of the "mean 
excitation energy" in {fi~ , e~) is different from that in z/^) | ]101| |, because 



in the last process the coherent channel doesn't exist. For the e~) process 
the coherent channel {gs gs) appears only in the denominator of eq. ( |199| ) 
and, since this dominates the total /i — e conversion rate, the resulting mean 
excitation energy E in this process is much smaller than that for the i/^) 
reaction. Obviously, the mean excitation energy E can be evaluated by 
constructing explicitly all the possible excited states | / > in the context of a 
nuclear model. In ref. ||101|| , for example, the QRPA method has been used 



for the determination of the mean excitation energy of the ^^Ti nucleus in 
the process /i~ Ti — »• e~ -|-^^ Ti*. 

The method of closure approximation proceeds by defining the operator taken 
as a tensor product from the single-particle operators T of eq. (|163|) or ( |164| ). 
For a O"*" initial (ground) state the relevant tensor product is 

O = J2[f'' X f*-'']° (200) 

k,k' ^ 

The corresponding total matrix elements in a sum-rule approach are written 

as 



ML 



r-^Y < I \ Ov \ t > + 3fi < t \ 6, 



I > 



(201) 



where | k | is the average momentum and the operators Oy (vector) and Oa 
(axial vector), which contain both one-body and two-body pieces, result from 
the corresponding T operators of eqs. (|163| ), ( |164^ ). Consequently, one has 
to evaluate the matrix elements of the operators Oy, Oa in a given model. 

3) Nuclear matter mapped into nuclei with local density approximation: 

The method has been already discussed in sects. 4.2.1-4.2.2. The obtained 
results will be presented in sect. 5.3. 
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5.2.1. RPA calculations envolving the final states explicitly. In ac- 
tual shell model calculations it is quite hard to construct the final states 
explicitly in medium and heavy nuclei. In such cases we can employ the RPA 
approximation. In the context of the quasi-particle RPA the final states en- 
tering the partial rate matrix elements are obtained by acting on the vacuum 



> with the phonon creation operator [|108i p. 101 , |114 



Q+(/JM) = 5^[Ay'^'"M+(a, JM) - Y^'-^'^^Aia, JM)] (202) 

a,T 

(angular momentum coupled representation) i.e. | / >= | >. The 
quantities X and Y in eq. ( P02| ) are the forward and backward scattering 
amplitudes. The index a, runs over proton (r = 1) or neutron (r = —1) two 
particle configurations coupled to J. 

The nuclear matrix element involved in the partial rate r^.^^ from an initial 
state 0"*" to an excited state | / > takes the form 



< / II f (''^)^ II 0+ >= ^ w^/ [ xy'^'^)[/j;V/;) + i-fy^'^'^^vj^^u^:^ 

(203) 

The phase ^=0, 1 for the vector, axial vector operator, respectively ||116|| . The 
probability amplitudes V and U for the single particle states to be occupied 
and unoccupied, respectively, are determined from the BCS equations and 
the X and Y matrices are provided by solving the QRPA equations. The 
quantities = W^^^^ contain the reduced matrix elements of the operator 
T between the single particle proton or neutron states ji and j2 as 

w-r ^<n\\f^\\n> ^ ^ 

In table 3, the results of ^^Ti are shown for various values of the parameter f3 
in a model space consisting of all the single particle levels included up to Shoo 
(same for protons and neutrons). In the photonic mechanism mechanism the 
axial vector matrix elements are zero (see eq. ( |160| )). 
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Table 3. Incoherent /i — e conversion matrix elements: vector {Sv), for the 
photonic mechanism (P = 3) and vector and axial vector {Sa), for a non- 
photonic mechanism ((3 = 5/6). They are for all the excited states included 
in the up to 3huj model space for ^^Ti. 



Mode photonic mechanism non-photonic mechanism 



Sv Sv Sa 



0+ 


1.111 


2.363 


0.0 


1+ 


0.0 


0.0 


0.297 


2+ 


0.309 


0.422 


0.046 


3+ 


0.0 


0.0 


0.050 


4+ 


0.002 


0.002 


2. 10-^ 


5+ 


0.0 


0.0 


2. 10-6 


6+ 


2. 10-6 


2. 10-6 


2. 10-^ 


0- 


0.0 


0.0 


0.818 


1- 


9.744 


17.853 


0.795 


2- 


0.0 


0.0 


0.670 


3- 


0.052 


0.068 


0.011 


4- 


0.0 


0.0 


0.010 


5- 


8. 10-5 


1. 10-^ 


1. 10-5 


6- 


0.0 


0.0 


1. 10-5 



Total 11.217 20.708 2.697 



5.2.2. Sum-Rules in the context of QRPA. The RPA sum-rules for 
the total (/i-,e-) rate, in the case when an uncorrelated ground state vac- 
uum I i >=| > in eq. (|201|) is used, can be easily obtained. Then the 
matrix elements < | O | > are given by 



< I O I > 



(205) 
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[Y^{2j + l)<j\\T^\\j>{v^ 



+ (2J+1)^ p(aJ,r)p(aJ,r) } 



where the first term gives the one-body contribution and the second the 
two-body contribution. The quantities p, p are given by 



p(a J, • 



mi 



(206) 



r{r)Tj(r) 
32 



(207) 



« = (j2 > ji)- The phase 6 is the same as in eq. (|203|) . 



The correlated quasi-particle RPA sum- rule of the (yU~, e~) reaction, which is 
the expectation value of the operator O of eq. ( pO0| ) in the correlated ground 
state I > of eq. ( |191| ), is written as 



< I 6 I > = iVo^ ( < I 6 I > + < I {[5, 6] + [6, 5+] } I > ) 

(208) 

where the second term in the brackets takes into account the 2p2h excitations. 
Thus, in order to compute the total rate with the correlated matrix elements, 

< I O I >, one should first calculate the uncorrelated matrix element 

< I O I > and afterwards the contribution coming from the ground state 
RPA correlations. The latter contribution can be cast in the form 



< 1 {[s,d] + [0,5+]} 1 > 



(209) 



J2]l J 2J 1 



ab,X.T 



;i + '^i2il)(l+ Wl) 



The summation (|210D gives the proton-proton (p-p) and neutron-neutron (n- 
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n) correlations. The results obtained for ^^Ti (see tables 4 and 5) show that 
the contribution from the two-body correllations as expected are small. 



Table 4- Total rate nuclear matrix elements and gs — > gs transition for the 
photonic mechanism (^^,e~) conversion rates in "^^Ti calculated: 1) with 
shell model sum-rule and 2) with various types of QRPA sum-rules for dif- 
ferent mean excitation energies. 



Method 






E 


ML 




Shell Model{sum — rule) 


144.6 


20.0 


188.8 


67.2 


QRPA {explicit) 




135.0 




161.0 


83.9 


QRPA{sum — rule) 




135.0 


1.7 


138.3 


97.6 


QRPA{sum — rule) 




135.0 


5.0 


140.6 


96.0 


QRPA{sum — rule) 




135.0 


20.0 


141.7 


95.3 


QRPA + Corr {sum 


— rule) 


87.8 


1.7 


90.4 


97.1 


QRPA + Corr {sum 


— rule) 


87.8 


5.0 


91.8 


95.6 


QRPA + Corr {sum 


— rule) 


87.8 


20.0 


92.6 


94.8 



Table 5. Total rate matrix elements and gs gs transitions for Ti given 
by various methods for the non-photonic mechanism P — 5/6 (see caption of 
table 4)- 



62 



Method 




M2 

Cj S — > Cf s 


E 


ML 

iOZ 


V (%) 


Shell Model{sum — rule) 


374.3 


20.0 


468.0 


80.0 


QRPA {explicit) 




363.0 




386.4 


93.9 


QRPA{sum — rule) 




363.0 


0.5 


366.2 


99.1 


QRPA{sum — rule) 




363.0 


5.0 


376.5 


96.4 


QRPA{sum — rule) 




363.0 


20.0 


382.8 


94.8 


QRPA + Corr {sum 


— rule) 


236.2 


0.5 


238.6 


99.0 


QRPA + Corr {sum 


— rule) 


236.2 


5.0 


245.4 


96.3 


QRPA + Corr {sum 


— rule) 


236.2 


20.0 


249.5 


94.7 



5.2.3 Sum-rules in the context of shell model. The total /i-e conver- 
sion rate matrix elements of eq. ( ^01| ) can be conveniently evaluated by using 
shell model closure approximation. To obtain the needed matrix elements we 
assume that the initial nuclear wave function | i > is a Slater determinant 
with closed proton and neutron (sub) shells constructed out of single particle 
harmonic oscillator wave functions. The tensor product operator O of eq. 
(|200|) in this case (no second quantization) is written as ||115| 



Oa = glT. + Ba{r3i + Ts,) + C^r^.T^j] a = V,A (210) 

where 



Ay = 9, Bv = 9, Cv =9, ~gy = \^~gA = Q (211) 

Av = 9//^, By = {3/fv)f3, Cv = /5^ gv = | 

Aa = e/fv, Ba = (e//y)/3, Ca = /5^ ~9a = | (212) 

for the photonic and non-photonic cases correspondigly. Furthermore, 

eS = jo(gr.,), e^.=jo(gr.,)^ (213) 



63 



From eq. ( |210| ) we can see that, contrary to the z/^) reaction, a distinct 
feature of the operator O, which is responsible for the process e~), is the 
presence of an isospin independent term in both the vector and axial vector 
component. It has been shown | |115| | that this part of the fi — e operator, is 
the dominant one. By using standard shell model techniques |3l| the matrix 



elements M^^^, which contain one body and two body terms, can be cast in 
compact analytical forms as 



S, = g,{A,Z){^) ^,a'^)e~'^'/\ a = | q | 6, X = V, A 

(214) 

where b is the harmonic oscillator parameter, N^ax is the maximum number 
of oscillator quanta occupied by the nucleons in the considered nucleus and 
are appropriate coefficients which may depend on A and Z. The functions 
gx{A, Z) describe the total rate if its dependence on the momentum transfer 
can be neglected. The definitions of the coefficients ^\ and the functions 



gx{A, Z), are given in ref. |3T 



The total (/i , e ) matrix elements using the shell model sum rules take the 
form 



Ml =~gl[z+ (|-^) V + A^ + ^(^) + Sj^\ 



a = V,A 



(215) 

where S"^, k = 0, 1, 2 correspond to the three operators entering eqs. ( piO| )- 
( 213 ). In tables 3 and 4 the results obtained this way in ^^Ti for average 



exitation energy E = 20MeV, are compared with those obtained by using 
quasi-particle RPA. 

5.2.4. The {fi~,e~) conversion in the PrimakofF's method. It is well 

known that success of the sum rule method hinges upon a good choise of 
mean excitation energy E. Then it requires only knowledge of the structure 
of the ground state. Very early Primakoff developed a phenomenological 



method for the (/x , t'^) reaction, which reproduces very well the experimental 
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data for the total muon capture rate and which does not contain the energy 



E. It yields the well-known Goulard-Primakoff ||6^ function (see sect. 3.2.2). 
which contains three parameters, the strengths of the isoscalar, isovector and 
isotensor parts of the muon-capture operator which are determined by a fit 
to the total muon capture reaction data of the hole periodic table. 

By exploiting the common components of the (/i~,z/^) and (/i~,e~) conver- 
sion operators one can construct a phenomeno logical formula for the e~) 
process [^. In this way the isospin dependent total fi — e rate can be esti- 
mated by using the values for the isoscalar, isovector and isotensor param- 
eters of the Goulard-Primakoff function determined by the ordinary muon 
capture data. Even though the results of this phenomenologilal method com- 
pare well with those of the shell model closure approximation, the fraction of 
the total (/i — > e) rate coming from the isospin dependent part of the e~) 
operator is tiny. Furthermore, as we have stressed in sect. 4, the dominant 
channel in the {fi~,e~) process is the coherent one which is not possible in 
the (/i, z/^) and therefore it cannot be obtained from the Goulard-Primakoff 
method. 



5.3. Discussion of nuclear matrix elements 

As it can be seen from table 4, the ground state to ground state transition 
(coherent contribution) exhausts a large portion of the sum rule given by the 
parameter rj. This is not entirely unexpected, since in this case the contri- 
bution of all nucleons is coherent. This is very encouraging since, as we have 
seen in sect. 2.3, the coherent process is of experimental interest. In order to 
extract r] in addition to the coherent rate we must have a reliable estimate of 
the total rate. The value of rj depends, of course, on the parameters Qv^Qa 
and (3 of the elementary amplitude. 

The only nuclear information needed for the coherent mode are the "muon- 
nuclear" form factors Fp{q^) and Fn{q^), (see eq. ( |175| )), which can easily 



be calculated both in the shell model and in RPA calculations. Admittedly, 
however, these calculations may not be very precise at the high momentum 
transfer = w?^c^ ~ 0.25/m~^. If on the other hand the muon wave function 
is assumed to be constant to be taken out of the radial integrals, the nuclear 
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matrix element essentially depends on the proton and neutron form factors 
F]s[{q^) and Fzip"^) which can be taken from experiment (electron scattering 
data, pionic atoms etc.). 



Using the results of table 6 (Chiang et ai, 1993 |99]), we can compute the 
parameter 



am,, . 



(216) 



This parameter gives us a measure of the deviation of the effective form factor 
(involving and the muon) from the usual nuclear form factor. The obtained 
results are presented in table 7. 



Table 7. Values of the parameter C,p.n{A, Z) of eq. in the text. The 

approximation of eq. ( 161 ) is not very accurate for heavy nuclei. 



(A,Z) 


(12,6) 


(32,16) 


(40,20) 


(48,22) 


(90,40) 


(208,82) 




0.979 


0.962 


0.946 


0.990 


0.850 


0.768 


in 


0.979 


0.956 


0.951 


0.639 


0.560 


0.285 



We thus see that, for heavy nuclei the muon wave function cannot be taken 
as a constant. This is especially true for the neutron component. This means 
that in this region the experimental form factors must be used with caution, 
even though the effect on the branching ratio may be less pronounced. 

All calculations of the coherent rate indicate that, in spite of the earlier 



expectations (Weinberg and Feinberg, 1959 [0, the coherent branching ratio 
increases all the way up to the Pb region where it starts decreasing. 
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For the supersymmetric model discussed in sect. 3.3.3, we present our results 
in table 8. 

Table 8. The nuclear form factors Fz and Fn entering in the quasi -elastic 
(/i, e) conversion (eq. ( \14^ ). The quantities •jph and k are defined in the 
text (see eqs. dlSOj ), l \151\ )), and are evaluated in the case of supersymmetric 
models. 



A 


Z 


Fz 


Fn 




K 




ReN 


/ Re^ 


4 


2 


0.865 


0.865 


1 


67 


1.56 


1.51 


10" 


-3 


12 


6 


0.763 


0.763 


1 


67 


3.64 


3.53 


10- 


-3 


14 


6 


0.753 


0.745 


1 


88 


7.96 


9.36 


10- 


-3 


16 


8 


0.736 


0.736 


1 


67 


4.52 


4.39 


10- 


-3 


28 


14 


0.639 


0.639 


1 


67 


5.95 


5.78 


10- 


-3 


32 


16 


0.618 


0.618 


1 


67 


6.37 


6.19 


10- 


-3 


40 


20 


0.582 


0.582 


1 


67 


7.05 


6.85 


10- 


-3 


48 


20 


0.563 


0.515 


1 


85 


16.08 


1.84 


10- 


-2 


48 


22 


0.543 


0.528 


1 


77 


9.74 


1.18 


10- 


-2 


60 


28 


0.489 


0.478 


1 


74 


9.24 


9.58 


10- 


-3 


72 


32 


0.456 


0.435 


1 


79 


11.54 


1.25 


10- 


^2 


82 


32 


0.440 


0.379 


1 


89 


24.98 


3.00 


10- 


-2 


88 


38 


0.412 


0.370 


1 


79 


12.98 


1.40 


10- 


-2 


90 


40 


0.406 


0.367 


1 


76 


11.41 


1.20 


10- 


-2 


114 


50 


0.335 


0.306 


1 


77 


10.35 


1.10 


10- 


-2 


132 


50 


0.315 


0.250 


1 


86 


25.80 


3.00 


10- 


-2 


156 


64 


0.263 


0.207 


1 


76 


11.96 


1.30 


10- 


-3 


162 


70 


0.253 


0.202 


1 


70 


8.92 


8.94 


10- 


'3 


168 


68 


0.249 


0.191 


1 


76 


12.47 


1.35 


10- 


-2 


176 


70 


0.242 


0.181 


1 


76 


13.75 


1.49 


10- 


-2 


194 


82 


0.198 


0.168 


1 


77 


7.20 


7.69 


10- 


-3 


208 


82 


0.189 


0.135 


1 


73 


10.42 


1.07 


10- 


-2 



The total rate can also be calculated in the context of the shell model. 
In this case it is not possible to construct all final states explicity in realistic 
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model spaces. One is thus forced to invoke the closure approximation (see 
sect. 5.2.3). One, however, has no idea about the proper average excitation 
energy Eexc to use. In the earlier calculations (Kosmas and Vergados, 1990 
E^) the value used was the same with that of the common muon-capture. 
Using this value of E^xc = 20MeV, the value of r] ranged from over 90% 
in light nuclei to about 30% in heavy nuclei (ibid. Table 4). We have seen, 
however, in sect. 5.2.2 that, with the proper definition of eq. ( |199| ), Eexc 
must be quite a bit lower. We should, therefore, give more credibihty to the 
RPA results discussed above {tables 4 and 5). 

The calculation of the total rates is quite a bit harder. It can, however, easily 
be done in the context of QRPA. In QRPA one can also apply the sum rule 
techniques with an average energy defined by eq. (|199|) , which in this case 
can be calculated. The agreement between the two methods is quite good 
(see tables 4 and 5). It is important to note that the ground state correlations 
tend to decrease all rates in our example ||101|| by 35%. The obtained value 
of r] is quite high. 

The total {fi, e) conversion rate has also been computed by a recent new 
method which utilizes nuclear matter mapped into nuclei with the local den- 
sity approximation (see sects. 4.2.1 and 4.2.2). This method has the advan- 
tage that both the incoherent rates of (/i, e) conversion and ordinary muon 
capture can be computed in the same way (see expressions for Gp{N, Z) and 
Gn{N, Z) of eqs. (p| ) and (p5|) as well as eq. (p^). Using the results of 



this calculation we see from table 6 that the parameter r] is in all cases greater 
than 80% and keeps increasing from light to heavy nuclei. Furthermore, the 
nuclear matrix elements for the coherent mode and the resulting branching 
ratios increase all the way to the Pb region. 



6. CONCLUDING REMARKS 

As we have seen in sect. 3, the most popular scenario for lepton-flavor vio- 
lation involves intermediate neutrinos at the one -loop levels. We have seen, 
however, that due to the GIM mechanism we encounter an unfavorable ex- 
plicit dependence on the neutrino mass. We have seen in sec. 3.2.1 that, 
for leptonic currents of the same helicity the amplitude for light neutrinos 
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depends on the square of the neutrino mass while for heavy neutrinos, on 
the inverse neutrino mass squared (see eqs. BO, PTI RBI, Bl). Thus in this 



mechanism lepton flavor is unobservable for neutrinos much lighter or much 
heavier than the W-boson mass. This unfortunately happens to be the case 
with most currently fashionable models. The explicit dependence on the 
neutrino mass is somewhat improved (linear for light neutrinos, inversely 
proportional for heavy neutrinos) in the case of L-R interference in the lep- 
tonic sector (see eqs. |8^, |8^ and |8^, ^). However, the situation essentially 
does not improve due to the presence of the mixings U^^'^^ or U^"^^^ which in 
most fashionable models is negligible. 

From the phenomenological point of view, therefore, we can add very little to 
the discussion of a previous review [12] and we will not elaborate further on 
this mechnism. We will instead summarize the results obtained using other 
intermediate particles. 

We will begin with an extended higgs sector involving non-exotic particles, 
i.e. two doublets in the Bjorken - Weinberg mechanism. The calculation 
of the branching ratio is rather complicated since the two loop contribution 
becomes dominant. It has recently been done by Chang, Hou and Keung 
[[70(1 , and Barr and Zee []7T| only for the /i — > 67 process. The predictions 
of this model depend essentially on four parameters: Au (the ttH Yukawa 
coupling), the combination Ae^coscpa, i^efj, is the efiH Yukawa coupling and 
cos(j)a (the W~^W~H cubic coupling relative to that of the stantard model), 
the top quark mass rrit and the Higgs scalar mass tuh- The authors take 
the somewhat optimistic choice of At* ~ 1 and Ae^cos0Q, = 1 and plot 
their obtained branching ratio as a function of and rriH- For ttih very 
large {niH 3> ITeV) or small {rriH <^ 200GeV) their results are essentially 
independent of rrit with branching ratios close to the present experimental 
limit. These authors, however, present their results as though there is one 
Higgs Boson by assuming that the contributions from the at least two needed 
Higgs Scalars do not cancel each other. 

Unfortunately, the above elaborate calculation does not provide the form 
factors Jeo and /mo which do not contribute to /i — 67 but enter into all the 
other lepton flavor violating processes which involve virtual photons. 
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The next extension of the standard model involves singly and/or doubly 
charged Higgs scalars. Lepton flavor violation can be caused by introduc- 
ing only the singly charged isosinglet both for fi ^ ej and fi —>■ 3e. The 
branching ratio unfortunately depends on a big power of the not accurately 
constrained mass of the isosinglet (the 4th for fi — > 67, the 8th for fi 3e). 
It also depends on the Yukawa couplings and Xer of the isosinglet. Thus, 
one cannot make accurate predictions. One instead can use the present ex- 
perimental limit for /i — >^ 67 to derive the constraint 



Ms > 94GeV{10\X,Aery^') (217) 

The introduction of the doubly charged Higgs scalars (isotriplet or isosinglet, 
see sect. 3) can give rise to lepton flavor violation at the tree level for fi ^ 3e 
and muonium-antimuonium oscillations which thus become favorable. In the 
case of the doubly charged isosinglet the amplitude for the /i 3e 
process can be written [|12] as 



^ Gp I- P1 2 
V2 V2 



M = n,^ ^ u{pe) (1 + l5Mp^)u{p2) (1 - I5MP1) (218) 



which leads to the branching ratio 

R=l\hx\^ (219) 

o 

with 

= (220) 

Once again the branching ratio depends on the inverse 4th power of the 
unknown mass of x^^. Assuming that g^fj, ~ ^fge — 0.1 and rrix = lO^GeV 
we obtain 

= 1.6 X 10~^°, i?~3.2xl0~^^ (221) 

which is many orders of magnitude away from the planned experiments. 
The situation with the doubly charged isosinglet is analogous except that 
75 —75 and fix ^ h^. Since, neither the mass not the couplings of the 
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isosiglet are determined from other experimental data, we can use the present 
experimental limit to constrain h/^. We obtain 



nA < 3 X 10"^ (222) 

which leads to the constrain 

^^4^ < 2 X 10-i°GeV^-' (223) 

The amplitude for muonium-antimuonium oscillations is the same with that 
of eq. ( pl8| ) except that, the antisymmetrization term (1 — P12) / a/2 is absent. 
Thus, for the isotriplet we see that the calculated value of is much smaller 
than the experimental limit < 0.16 obtained in sect. 2.4. Similarly, the 
limit fi-A < 3 X 10~^ extracted from the experimental limit of yU — > 3e is 
much smaller than 0.16. This, of course, indicates that if the doubly charged 
Higgs scalars exist it is much more likely to observe lepton flavor violation 
in /i — > 3e rather in M — M oscillations. 

We will finally discuss lepton flavor violation in supersymmetric theories (see 
sects. 3.3.1 and 3.3.3). For /i — > 67 the relevant equations are ( |134| ) and 
( |135| ). Assuming that m^i2 = — l50GeV, we obtain 



Re-y-8x IQ-^^ (224) 
Similarly, for the n ^ 3e the branching ratio we obtain 

i?3e ~ 5 X 10"^^ (225) 
i.e. /i — >^ 67 is favorable in this model. 

Let us now discuss (/i, e) conversion in supersymmetric theories. In addition 
to the parameter fj of (|134|) , which contains both the effect of renormalization 
and the mixing angles, we encounter all three functions f{x), g{x) and fb{x) 
we met in /x — 3e. For the experimentally interesting coherent contribution 
the branching ratio is given by eq. (|149|) . By noting that = — and 
rriu ~ rrie ~ "013/2 and assuming further that the photino is the lightest 
particle, one is led to eq. ( |154| ). In summary then one can write 
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ReN CX ,1 



+ -;f^)\h (226) 



All the relevant nuclear physics is contained in the parameters k and ^yph- 
These parameters are given in table 7. For comparison purposes we mention 
that for the nucleus *^^Ti, which is of experimental interest, the form factors 
become Fz = 0.538, Fn = 0.506 (QRPA) and Fz = 0.550 and F^ = 0.520 
(Fermi distribution). For the reader's convenience we also present in table 7 
the ratio ReN/Re-y- From this table we see that the coherent effect of all nu- 
cleons, tends to enhance the ratio R^n / Re-y compared to the naive estimates 
mentioned in the introduction. Wc see that, ReN still remains about two 
orders of magnitude less that R^j. This, however, can be compensated by 
the desirable experimental signature of (/x, e) conversion mentioned in sec. 
2.3. 

It is apparent that lepton flavor, unlike strangeness, may be absolutely con- 
served if the Gods of physics bestowed upon the standard model of elec- 
troweak interactions absolute authority. Very few people, however, subscribe 
to this dogma. Lepton flavor violation follows naturally in most extensions 
of the standard model. Its observation, though, is not going to come easy. 
The experimental efforts have reached limits which make further improve- 
ments extremely difficult. Worst yet the predictions of currently fashionable 
models are not encouraging them. So, barring unforseen developments, such 
efforts may be classified in the pursuit of nothingness. This, under the dif- 
ficult conditions of present economies may be catastrophic. But the key hes 
in "unforseen circumstances". Historically, this has been not only the rule 
but the main buty of science. In any case the theoretical predictions refiect 
present biases and should not deter the continuation of experimental efforts*. 



i/ uviy TO iirjSeu, Km acev o Koa/ioq o /iLKpoq o /le^aq. (227) 
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